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Abstract 

Quantum channels, also called quantum operations, are linear, trace 
preserving and completely positive transformations in the space of quan- 
tum states. Such operations describe discrete time evolution of an open 
quantum system interacting with an environment. The thesis contains an 
analysis of properties of quantum channels and different entropies used to 
quantify the decoherence introduced into the system by a given operation. 

Part U of the thesis provides a general introduction to the subject. In 
Part [U the action of a quantum channel is treated as a process of prepa- 
ration of a quantum ensemble. The Holevo information associated with 
this ensemble is shown to be bounded by the entropy exchanged during 
the preparation process between the initial state and the environment. A 
relation between the Holevo information and the entropy of an auxiliary 
matrix consisting of square root fidelities between the elements of the en- 
semble is proved in some special cases. Weaker bounds on the Holevo 
information are also established. 

The entropy of a channel, also called the map entropy, is defined as 
the entropy of the state corresponding to the channel by the Jamiolkowski 
isomorphism. In Part II I II of the thesis, the additivity of the entropy of 
a channel is proved. The minimal output entropy, which is difficult to 
compute, is estimated by an entropy of a channel which is much easier to 
obtain. A class of quantum channels is specified, for which additivity of 
channel capacity is conjectured. 

The last part of the thesis contains characterization of Davies channels, 
which correspond to an interaction of a state with a thermal reservoir in 
the week coupling limit, under the condition of quantum detailed balance 
and independence of rotational and dissipative evolutions. The Davies 
channels are characterized for one-qubit and one-qutrit systems. 
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Part I 

Introduction 

1 Preliminary information 
1.1 Preface 

It is not easy to give a satisfactory definition of information in sense in 
which this word is used in everyday life. For instance one could ask, 
how much information is contained in an allegorical baroque painting of 
Vermeer. There exist, of course, many interpretations and therefore, many 
kinds of information concerning this picture. However, nowadays we are 
willing to distinguish some sort of information necessary to communicate a 
message independently on the interpretation. Due to our experience with 
computers we are used to problems how to encode the information into 
a string of digital symbols, transmit it and decode it in order to obtain 
the original message in another place. Imagine that we need to send the 
information contained in the Vermeer's picture. We have to encode it 
into digital data, transmit it to the other place and recover the picture on 
the screen of the receiver's computer. In a sense we send almost all the 
information without knowing what interpretations it may carry. 

The problem rises what is the minimal amount of information mea- 
sured in binary digits that enable the receiver to reliably recover the orig- 
inal message. In considered example we can divide the image of the Ver- 
meer's picture into small pieces, decode colours of each piece into digital 
strings and transmit the description of colours one after another. How- 
ever, we can also save some amount of digits when we menage to describe 
shapes of regions of the same colours in the picture and send only infor- 
mation about colours, shapes and patterns. How to do that in the most 
efficient way? This is a major problem for experts working on the infor- 
mation theory and computer graphics. Some rules of the optimal coding 
were used intuitively during construction of the Morse alphabet. The let- 
ters which occur in the English language more frequently are encoded by 
a smaller amount of symbols. 

In communication and computer sciences the problem of data compres- 
sion is a subject of a great importance. To what extend the data can be 
compressed to still remain useful? Claude Shannon worked on the prob- 
lem of transmission of messages through telecommunication channels. In 
1958 he published his famous paper II opening the new branch of knowl- 
edge known as the theory of information. In this theory a message is 
composed of letters occurring with specified frequencies related to prob- 
abilities. Every letter of a message can be encoded as a string of digital 
units. Every digital unit can appear in one of r possible configurations. 
Shannon found what is the minimal average amount of digital units per 
symbol which encodes a given message. This smallest average number of 
digital units is related to the information contained in the message and is 
characterized by a function of the probability distribution P = {pi, ...,p„} 
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of letters, now called the Shannon entropy, 

n 

ff(p) = -X> lo g,Pi, (i) 

1 = 1 

where 01og r = 0, n is a number of letters and the base of the logarithm 
r characterizing the amount of configurations of a chosen digital unit can 
be chosen arbitrary. If the base is equal to 2, the unit of entropy is called 
binary unit or bit. 

The idea of efficient coding concerns in replacing more frequent letters 
by means of a smaller amount of bits. Shannon treated the message as a 
sequence of letters generated independently according to the probability 
distribution P specified for a given language. The original reasoning of 
Shannon proceeds as follows. There are so many possible messages as the 
amount of typical sequences of letters with a given probability distribution 
in the string of length k — > oo. Atypical sequences such as strings of letters 
a repeated k times are unlikely and are not taken into account. The 
amount of possible messages is given by the amount of typical sequences, 
which is of order of 2 kH ( p ^ if the base of the logarithm is equal to 2. This 
number is justified by methods of combinatorics. Hence, every typical 
message of length k can be represented by a string of bits of size kH(P). 
Therefore, the entropy H(P) can be interpreted as the smallest average 
amount of bits per letter needed to reliably encode each typical message. 

The information theory treats, as well, the information as a measure 
of uncertainty about the outcome of a random experiment. Looking for a 
function which is suitable as a measure of the uncertainty about the con- 
crete result of experiment, provided the probabilities of all experimental 
outcomes are given, Shannon formulated a few postulates for the infor- 
mation measure [T|: 

• It is a continues function of the probability distribution. 

• If all events are equally likely the function of uncertainty is an in- 
creasing function of their number. 

• If one of the events is split into two, the new function of uncertainty 
is equal to the sum of the original uncertainty and the uncertainty 
of the new division weighted by the probability of the divided event. 

The only function which satisfies these postulates is the Shannon entropy 
H(P). Therefore, the uncertainty or lack of information on the outcome 
of an experiment is the second interpretation of the entropy H(P). 

Taking a weaker set of axioms allows one to generalize the definition 
of the measure of uncertainty and to find other functions of probability 
vector P, which in special case converge to the Shannon entropy Q. For 
instance, Renyi introduced one parameter family of generalized entropy 
functions. Since, the information of an experiment consisting of two inde- 
pendent experiments should be given by the sum of the information gained 
in both experiments, the measure of information should be additive. The 
Shannon entropy of the joint probability distribution of two independent 
variables is additive. Renyi noticed [2] that the additivity of information is 
not equivalent to the third postulate of Shannon. However, if one replaces 
the third postulate by additivity of information of independent events, yet 
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another axiom should be postulated to obtain back the Shannon's formula 
This additional postulate specifies the way of calculating the mean 
values. If one considers the linear mean, the Shannon entropy is singled 
out by this set of postulates. However, other definition of the mean value 
also can be taken. In consequence, the new set of postulates implies an 
one parameter family of generalized entropy functions known as the Renyi 
entropy of order q: 



Here, q denotes the free parameter depending on the definition of the av- 
erage. Another generalization of entropy function was analysed by Tsal- 
lis [3]|4]. The Tsallis entropy of order q is defined as follows, 



Hence the information theory concerns entropies, however, it also in- 
vestigates communication sources and communication channels which can 
introduce some errors to messages. Information theory defines such quan- 
tities as the relative entropy and the mutual information [T]. Using these 
concepts the channel capacity is defined. It is the maximal rate of infor- 
mation which can be reliably decoded after passing through the channel. 
The channel capacity is measured in bits per a unit of time. 

The theory of quantum information, which considers quantum systems 
as carriers of information, should enable one to generalize the notions of 
classical information theory such as the channel capacity. To formulate 
a quantum counterpart of the Shannon concepts such as the relative en- 
tropy or channel capacity the theory of open quantum systems, quantum 
statistical processes, statistical operators, density matrices, partial traces 
and generalized measurements should be applied. In the early stage of the 
stochastic theory of open quantum systems, it was developed by Davies [5J, 
and Kossakowski [Bj. Moreover, other important results on accessible in- 
formation transmitted through a noisy quantum channel were obtained 
by Holevo [7|. 

There are many advantages of using quantum resources to transform 
and transmit the information 8[. In particular, there exist a famous 
protocols of superdense coding [9] of information into a quantum carrier. 
Furthermore, some computational problems can be solved in framework of 
the quantum information processing faster than classically |10H12| . Quan- 
tum world gives also new communication protocols like quantum telepor- 
tation [9"lll3| which is possible due to quantum entanglement [14, 15] . In 
quantum case, entangled states can increase the joint capacity of two chan- 
nels with respect to the sum of the two capacities |16H18| . Also a new 
branch of cryptography was developed due to the quantum theory [19] . 
Although, these new possibilities are promising, manipulation of quan- 
tum resources is difficult in practice. In particular, the quantum states 
carrying the information are very sensible to noise, which can completely 
destroy quantum information. Moreover, probabilistic nature of quantum 
theory does not allow us to extract uniquely the information from quan- 
tum sources. Many restrictions and laws of quantum information theory 






(3) 
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are formulated in terms of inequalities of quantum entropies. The most 
significant quantum entropy is the one defined by von Neumann [201121] . 
which is a counterpart of the Shannon entropy. However, the other quan- 
tum entropies such like the Renyi entropy [2] or Tsallis entropy are also 
considered [51131122]. 

The issue of transmitting a classical information through a noisy chan- 
nel is an important issue in the information theory [1,23, 24] . Among 
problems concerning the information channels one can specify the fol- 
lowing questions: How to encode the information in order to transmit it 
reliably through the channel in the most efficient way 1,25]? What is the 
maximal rate of the information transmission? What is the capacity of a 
given communication channel 26-29J? Which states are the most resis- 
tant to errors occurring in the a channel [30,311? What are the efficient 
strategies of the error correction [32]? 

Similar questions can also be formulated in the framework of quantum 
information theory. The quantum channels, also called quantum oper- 
ations, are transformations in the set of states [33H36] . They describe 
evolution of an open quantum system interacting with an environment in 
discrete time. 

The set of all quantum channels is still not completely understood. 
Merely the set of one-qubit channels is satisfactory explored [37138] . How- 
ever, even in this simplest case some interesting problems are open. For 
instance, it is not clear, whether the capacity of one-qubit channels is 
additive |18| . Another approach to quantum channels suggests to analyse 
only certain special classes of them, motivated by some physical mod- 
els [3HH13]- 

Quantum channels are also formally related to measurement processes 
in quantum theory [351145] . As a measurement process changes the quan- 
tum state and in general cannot perfectly distinguish measured states, 
there is a fundamental restriction on the information which can be ob- 
tained from the message encoded into quantum states [7|. These restric- 
tions are also formulated in terms of entropies. 

The different aspects of quantum channels mentioned above suggest 
that entropies which characterize the channels play an important role in 
the information theory. This thesis is devoted to investigation of quan- 
tum channels and some entropies used to characterize them: the minimal 
output entropy [181139] . the map entropy 46-48] and the exchange en- 
tropy [25]. 

1.2 Structure of the thesis 

The thesis is partially based on results already published in articles [461491 - 
153] . which are enclosed at the end of the thesis. In a few cases some issues 
from these papers are discussed here in a more detailed manner. The thesis 
contains also some new, unpublished results and technical considerations 
not included in the published articles. 

The structure of the thesis is the following. The thesis is divided into 
three parts. The first part is mostly introductory and contains a short 
review of the literature. This part provides basic information useful in 
the other parts of the thesis and fixes notation used in the entire work. In 



part U only the result from Section \l . 6. 1 1 concerning the Kraus representa- 
tion of a complementary channels and Section fO|] on the Kraus operators 
constructed for an ensemble of states are obtained by the author. 

Part [IT] contains results based on papers [46, 49 , 521, not known before 
in the literature. However, some results not published previously are also 
analysed there. 

Chapter [2] contains the most important result of the thesis - the in- 
equality between the Holevo information related to an ensemble of quan- 
tum states and the entropy of the state of environment taking part in 
preparation of the ensemble. As the entropy of the environment can be 
treated equivalently as the entropy of an output of the complementary 
channel, or the entropy of a correlation matrix, or the entropy of a Gram 
matrix of purifications of mixed states, or as the entropy exchange, this 
relation might be considered as a new and universal result in the theory 
of quantum information. Consequences of this inequality have not been 
analysed so far. Chapter [2] contains also the discussion of the particular 
cases for which the inequality is saturated. This result has not been pub- 
lished before. Section 12.11 describes proofs of known entropic inequalities 
which are related to the bound on the Holevo quantity. Some new and 
unpublished consequences of these inequalities are presented in Section 
12.1.11 Original, new results are also contained in Sections 12.21 and 12.31 

Part [Tl] contains, moreover, the conjecture on the inequality between 
the Holevo information of a quantum ensemble and the entropy of the 
matrix of square root of fidelities. Several weaker inequalities are analysed 
here in a greater detail than it was done in [52] . Section [3.21 presents a 
confirmation of the conjecture for a special class of ensembles of quantum 
states. 

Part IIHI of the thesis is based on the results presented in [501151] . 
Article |51| described partially in Chapter|3]considers the relation between 
minimal output entropy and the map entropy. Section [4~2l contains a proof 
of additivity of the map entropy with respect to the tensor product of two 
maps, already published in our work |51| . These results allow us to specify 
a class of quantum channels for which additivity of the minimal output 
entropy is conjectured. 

The Davies maps acting on one-qubit and one-qutrit quantum systems 
are analysed in Chapter [5] Conditions for the matrix entries of a quantum 
operation representing a Davies map are given along the lines formulated 
in our work [50 j . Multiplicativity of the maximal output norm of one- 
qubit Davies maps, entirely based on the analogical proof for bistochastic 
maps 54 , is presented in Section 15.61 However, this result cannot be 
treated as a new one, since multiplicativity of the maximal output two 
norm was proved earlier for all one-qubit quantum channels [18| . Section 
l5.7l contains graphical representations of stochastic matrices of order three 
which correspond to quantum Davies maps, which has not been published 
yet. 

1.3 A short introduction to quantum mechanics 

The formalism of quantum mechanics can be derived from a few postulates 
(axioms) which are justified by experiments. The set of axioms defining 
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the quantum theory differs depending on the author 55J. However, some 
features occur common in every formulation, either as axioms or as their 
consequences. One of such key features is the superposition principle. It 
is justified by several experimental data as interference pattern in double 
slit experiment with electrons or interference of a single photon in the 
Mach-Zender interferometer |56| . The superposition principle states that 
the state of a quantum system, which is denoted in Dirac notation by \ip), 
can be represented by a coherent combination of several states \tpi) with 
complex coefficients a,i, 

|V>>=X>hk>. (4) 

i 

The quantum state \ip) of an N level system is represented by a vector 
from the complex Hilbert space Hn- The inner product {ipi\ip) defines 
the coefficients a* in Q. The square norm of ai is interpreted as the 
probability that the system described by \ip) is in the state \tpi). To 
provide a proper probabilistic interpretation a vector used in quantum 
mechanics is normalized by the condition {ip\ip) = ||^/>|| 2 = l a *| 2 — !• 

Quantum mechanics is a probabilistic theory. One single measurement 
does not provide much information on the prepared system. However, 
several measurements on identically prepared quantum systems allow one 
to characterize the quantum state. 

A physical quantity is represented by a linear operator called an ob- 
servable. An observable A is a Hermitian operator, A — A\ which can 
be constructed by a set of real numbers Xi (allowed values of the phys- 
ical quantity) and a set of states determined by the measurement, 
A — Y^,i Xi\ifii){ipi\. The physical value corresponds to the average of the 
observable in the state \ip), 

W^E^WI^W). (5) 

i 

One can consider the situation in which a state is not known exactly. 
Only a statistical mixture of several quantum states \<j>i) which occur with 
probabilities pi is given. In this case the average value of an observable 
has the form 

(A) {Pi , M =^>(&Wi>, (6) 

i 

which can be written in terms of an operator on Hn called a density 
matrix p = J2iPi\ < Pi){4>i\ as 

Ww,}=T*M- ( 7 ) 

A density matrix describes a so called mixed state. In a specific basis 
the density matrices characterizing an N level quantum system are rep- 
resented by N x TV matrices p which are Hermitian, have trace equal to 
unity and are positive. Let us denote the set of all such matrices by Mn, 

Mn = {p ■ dimp = N,p = p\p > 0,Trp = 1}. (8) 

This set is convex. Extremal points of this set are formed by projectors 
of the form \tf>) (if) | called pure states, which correspond to vectors \tp) of 
the Hilbert space. 
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The state of composed quantum system which consists of one iVi-level 
system and one A^-level system is represented by a vector of size N1N2 
from the Hilbert space which has a tensor product structure, Hn ± N2 = 
Hnx ®Hn 2 - Such a space contains also states which cannot be written 
as tensor products of vectors from separate spaces, 

IV12) * hM ® \i>2). (9) 

and are called entangled states. States with a tensor product structure are 
called product states. If the state of only one subsystem is considered one 
has to take an average over the second subsystem. Such an operation is 
realized by taking the partial trace over the second subsystem and leads 
to a reduced density matrix, 

pi = Tr 2 p 12 . (10) 

A density matrix describes therefore the state of an open quantum system. 

The evolution of a normalized vector in the Hilbert space is determined 
by a unitary operator \ip') — U\tp). The transformation U is related to 
Hamiltonian evolution due to the Schrodinger equation, 

ih±\i,)=H\i>), (11) 

where H denotes the Hamiltonian operator of the system, while t repre- 
sents time and 2nh~ is the Planck constant. A discrete time evolution of an 
open quantum system characterized by a density operator p is described 
by a quantum operation which will be considered in Chapter 11.61 

According to a general approach to quantum measurement [351157] . it 
can be defined by a set of k operators {_E l }* =1 forming a positive operator 
valued measure (POVM). The index i is related to a possible measurement 
result, for instance the value of the measured quantity. The operators E l 
are positive and satisfy the identity resolution, 

k 

53^ = 1 (12) 

i=l 

The quantum state is changing during the measurement process. After the 
measurement process that gives the outcome i as a result, the quantum 
state p is transformed into 

Pi = K l pK lt I Tr{K l pK lt ), (13) 

where K' — E % > 0. The probability p, of the outcome i is given by 
Pi = Tr(K * pK *'). Due to relation (|12|) . the probabilities of all outcomes 
sum up to unity. 

A quantum state characterizing a 2-level system is called qubit and its 
properties are discussed in more detail in Section fl. 71 

1.4 Schmidt decomposition 

The theorem known as Schmidt decomposition [58] provides a useful rep- 
resentation of a pure state of a bi-partite quantum system. 
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Theorem 1 (Schmidt). Any quantum state ^12} from the Hilbert space 
composed of the tensor product of two Hilbert spaces Hi®H2 of dimensions 
di and di, respectively, can be represented as 



|^12>=X)A i |tl>®|t2>, (14) 
i=l 

where { | * 1 ) } ^-L 1 and {|i2)}f=i are orthogonal basis of the Hilbert spaces Hi 
and H2 respectively, and d = min{di, d2}. 

Proof. Choose any orthogonal basis {|0i)}^Lj of Hi and any orthogonal 
basis {|02)}j=i °f ^2- In this product basis, the bi-partite state ^12) 
reads 

|Vl2)= Y, «wl*i>®l#>- (15) 

0<fc<di, 0<j<d 2 

Singular value decomposition of a matrix A of size di x d2 with entries 
dkj gives a k j — ^2 i u ki \ i Vij. Here u ki and Vij are entries of two unitary 
matrices, while Aj are singular values of A. Summation over indexes k 
and j cause changes of two orthogonal bases into 

|ii> (16) 

k 

\i2)=Y^Vij\4). (17) 

j 

The number o nonzero singular values is not larger than the smaller one 
of the numbers (di , CZ2) . □ 

The Schmidt decomposition implies that both partial traces of any 
bi-partite pure state have the same nonzero part of the spectrum: 

d 

Tr 1 \^ 12 ){^ 12 \ = ^2\ 2 i \i 2 ){i 2 \, (18) 
i=i 

d 

T* 2 |iM(iM = X)A?|ti><ii|. ( 19 ) 

i=l 

The Schmidt coefficients A, are invariant under local unitary transforma- 
tions Ui <8> U 2 applied to \tpi 2 ). The number of non-zero coefficients A, is 
called the Schmidt number. Any pure state which has the Schmidt num- 
ber greater than 1 is called entangled state. A pure state for which all 
Schmidt coefficients Aj are equal to 1/ \fd is called a maximally entangled 
state. 

Another important consequence of the Schmidt decomposition is that 
for any mixed state p there is a pure state \ip) of a higher dimensional 
Hilbert space such that p can be obtained by taking the partial trace, 

p = ThM<# (20) 

Such a state |V>) is called a purification of p. The Schmidt decomposition 
gives the recipe for the purification procedure. It is enough to take square 
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roots of eigenvalues of p in place of Xi and its eigenvectors in place of 
Any orthogonal basis in H2 provides a purification of p, which can 
be written as 

|^>=X)(^i®Vfl)|ii>®|ia>, (21) 

i 

where Uj is an arbitrary unitary transformation and v/pI* 2 ) — A»|ia). 

1.5 Von Neumann entropy and its properties 

Many theorems concerning the theory of quantum information can be 
formulated in terms of the von Neumann entropy [59 of a quantum state, 

S{p) = -Trplogp, (22) 

which is equivalent to the Shannon entropy |T|) of the spectrum of p. The 
entropy characterizes the degree of mixing of a quantum state. Assume 
that p is a density matrix of size N. The value of S(p) is equal to zero if 
and only if the state p is pure. It gains its maximal value log N for the 
maximally mixed state p* = il only. 

Von Neumann entropy has also an important interpretation in quan- 
tum information theory, as it plays the role similar to the Shannon entropy 
in the classical theory of optimal compression of a message |25| . Let the 
letters i of the message, which occur with probabilities pi, be encoded 
into pure quantum states from the Hilbert space Hn- Sequences 
of k letters are encoded into a Hilbert space of dimension N k . A long 
message can be divided into sequences of size k — > 00. Among them one 
can distinguish sequences in typical subspaces and such which occur with 
negligible probability. A unitary transformation applied to the sequence 
of quantum systems can transmit almost all the information into a typical 
subspace. The space of a typical sequence has the smallest dimensionality 
required to encode the message reliably with negligible probability of an 
error. This smallest dimensionality per symbol is shown [25| to be equal 
to the von Neumann entropy of the state p — Therefore, 
quantum coding consists in taking states from the smaller subspace of 
dimension 2 fcS ' p ' instead of a space of dimension N k to encode the same 
message. If the state p represents completely random set of states there 
is no possibility to compress the message, since S(p) = S(p*) = log 2 N, 
where logarithm is of base 2. The entropy, therefore, describes the capa- 
bility of compression of the message encoded in a given set of states, or 
the smallest amount of qubits needed to transmit a given message. 

The von Neumann entropy, as the entropy of eigenvalues of a density 
matrix, describes also the uncertainty of measuring a specific state from 
the set of the eigenvectors. The most important properties of the von 
Neumann entropy are |20| : 

• The von Neumann entropy is a non negative function of any p. 

• It is invariant under unitary transformations, S(p) = S(UpU^). 

• It is a concave function of its argument, Y^i=i ViS{pi) < S(%2i=i PiPOi 
where Pi > for any i and J^. pi = 1. 
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• It is subadditive 

S(jH2) < S( P i) + S(p 2 ), (23) 
where pi2 is a bi-partite state of a composite system and the partial 
traces read pi = TY2 pi% and p2 = Tri pi2. 

• The von Neumann entropy satisfies the relation of strong subaddi- 
tivity [BD], 

S(pi 23 ) + S( P2 ) < S(fiia) + S{p 23 ), (24) 
where the state P123 is a composite state of three subsystems (1, 2, 3) 
and the other states are obtained by its partial traces. 

1.6 Quantum channels and their representations 

One distinguishes two approaches to describe time evolution of an open 
quantum system. One of them starts from a concrete physical model 
defined by a given Hamiltonian which determines the Shrodinger equation 
(|11[) or the master equation, |45| . Solving them one may find the state of 
the quantum system at any moment at time. An alternative approach to 
the dynamics of an open quantum system relies on a stroboscopic picture 
and a discrete time evolution. It starts from a mathematical construction 
of a quantum map, p' = &(p), allowed by the general laws of quantum 
mechanics. This approach is often used in cases in which the physical 
model of the time evolution is unknown. This fact justifies the name 
"black box" model to describe the evolution characterized by a quantum 
map <E>. Such a model is also considered if one wants to investigate the set 
of all possible operations independently on whether the physical context 
is specified. Main features and some representations of the map <!>, which 
describes a "black box" model of non-unitary quantum evolution, are 
given below. 

The quantum map $ describes the dynamics of a quantum system p 
which interacts with an environment. It is given by a nonunitary quan- 
tum map $ : p — > p' . Any such map is completely positive, and trace 
preserving |33H36| . "Complete positivity" means that an extended map 
$ <S> 1m, which is a trivial extension of $ on the space of any dimension 
M, transforms the set of positive operators into itself. A completely pos- 
itive and trace preserving quantum map is called quantum operation or 
quantum channel. 

Due to the theorem of Jamiolkowski |34| and Choi |33| the complete 
positivity of a map is equivalent to positivity of a state corresponding to 
the map by the Jamiolkowski isomorphism. This isomorphism determines 
the correspondence between a quantum operation <E> acting on N dimen- 
sional matrices and density matrix D$ /N of dimension iV 2 which is called 
Choi matrix or the Jamiolkowski state 

Id* = [idiv®<E-](|0 + ><<^ + |), (25) 

where \<f> + ) = -7= 53 <=1 \i ) ® \i ) is a maximally entangled state. The dy- 
namical matrix D<s> corresponding to a trace preserving operation satisfies 
the partial trace condition 

Tr 2 -D$ = 1 (26) 
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The quantum operation $ can be represented as super-operator ma- 
trix. It is a matrix which acts on the vector of length N 2 , which contains 
the entries pij of the density matrix ordered lexicographically. Thus the 
superoperator $ is represented by a square matrix of size iV 2 . The super- 
operator in some orthogonal product basis ® is represented by a 
matrix indexed by four indexes, 

= (A ® 01*l*> ® 10- (27) 

The matrix representation of the dynamical matrix is related to the su- 
peroperator matrix by the reshuffling formula [15] as follows 

{i\ <S> (j\D*\k) <S> 10 = (i\ ® (k\<f>\j) ® 10- ( 28 ) 

To describe a quantum operation, one may use the Stinespring's dila- 
tion theorem [61 j . Consider a quantum system, described by the state p 
on Hn, interacting with its environment characterized by a state on Hm- 
The joint evolution of the two states is described by a unitary operation 
U. Usually it is assumed that the joint state of the system and the en- 
vironment is initially not entangled. Moreover, due to the possibility to 
purification the environment, its initial state is given by a pure one. The 
evolving joint state is therefore: 

ui = u(\l){l\®p)rf, (29) 

where |1) £ Hm and U is a unitary matrix of size JVM. The state of the 
system after the operation is obtained by tracing out the environment, 

M 

p' = $(p) = Tr M [ul\l}{l\®p)U t ] =J2 K *P K ' f > ( 3 °) 

i = l 

where the Kraus operators read, K l — {i\ U |1). In matrix representation 
the Kraus operators are formed by successive blocks of the first block- 
column of the unitary evolution matrix U. Here the state oj can be equiv- 
alently given as 

M 

u, = Ri p Rji ® io oi ■ (3i) 

A transformation p — > uj is obtained by an isometry F : Hn — > Hnm, 
where 

F|^=£(iT|0»®|;>. (32) 

i 

Due to the Kraus theorem |35| any completely positive map $ can be 
written in the Kraus form, 

AI 

1 = 1 

The opposite relation is also true, any map of the Kraus form (|33|) is 
completely positive. 



15 



1.6.1 Representation of a complementary channel 

Consider a quantum channel $ described by the Kraus operators K z , 

M 

Hp) = Tr M w = J2 Ki P Kit > ( 34 ) 

where notation from Section 11.61 is used. The channel <3> complementary 
to $ is defined by 

JV 

4(p) = Trjvo; = ^ K* pK^ (35) 

i=l 

and it describes the state of the M-dimensional environment after the in- 
teraction with the principal system p. One can derive the relation between 
operators {K-'}f =1 and {K'}ft 1 from the last equation by substituting ui 
as in (I31|) . This relation can be rewritten as 

AI JV 

(Tr K l pK^) \i) 01 = E (36) 
Comparison of the matrix elements of both sides gives 

N JV 
^ ^ Kj rnj pmnK n j = ^ ^ KgmPmnKna: (37) 

where matrix elements are indicated by lower indexes and the Einstein 
summation convention is applied. Hence, for any quantum channel $ 
given by a set of Kraus operators K l , one can define the Kraus operators 
K a representing the complementary channel <E> as 

K?j = i = l,...,M, j,a = l,...,N. (38) 



1.7 One— qubit channels 

One-qubit channels acting on density matrices of size 2 have many special 
features which cause that the set of these channels is well understood 
37,38,54 . However, many properties of one-qubit maps are not shared 
with the quantum maps acting on higher dimensional systems. Since one- 
qubit quantum channels are often considered in this thesis, the following 
section presents a brief review of their basic properties. 

A quantum two level state is called quantum bit or qubit. It is repre- 
sented by a 2 x 2 density matrix. Any Hermitian matrix of size two can 
be represented in the basis of identity matrix and the three Pauli matrices 

a = {(71, <T2, 0-3}, 

<n=(J J), «*.= (J J ), <r 3 =Q (39) 

One qubit state p decomposed in the mentioned basis is given by the 
formula 

1 -i 
p= -(id+f -a), few. (40) 
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Positivity condition, p > 0, implies that |f| < 1. The vector r is called 
the Block vector. All possible Bloch vectors representing quantum states 
form the Block ball. Pure one-qubit states form a sphere of radius |f| = 1. 

Any linear one-qubit quantum operation $ transforms the Bloch ball 
into the ball or into an ellipsoid inside the ball. The channel $ transforms 
the Bloch vector r representing the state p into r ' which corresponds to 
p' . This transformation is described by 

r' = Wr + K. (41) 

Here the matrix W is a square real matrix of size 3. A procedure analogous 
to the singular value decomposition of the matrix W gives W = O1DO2, 
where Oi represents an orthogonal rotation and D is diagonal. Up to two 
orthogonal rotations, one before the transformation $ and one after it, 
the one-qubit map $ can be represented by the following matrix 



/ 1 








°\ 


«1 


m 













<n 





\K3 








V3/ 



(42) 



The absolute values of the parameters rji are interpreted as the lengths of 
the axes of the ellipsoid which is the image of the Bloch ball transformed 
by the map. The parameters K4 form the vector ft of translation of the 
center of the ellipsoid with respect to the center of the Bloch ball. 

Due to complete positivity of the map $ and the trace preserving 
property, the vectors rf and k are subjected to several constraints. They 
can be derived from the positivity condition of a dynamical matrix given 
by HUEZ]: 



/ 1 + m + k 3 


«1 — i * K2 

\ rji + 772 





1 - 773 + «3 
T7l - 7?2 

tl — i * K2 



K\ + i * At2 
1 — 7)3 — «3 





771 + m 

Kl + i * K2 


1 + 773 - K3 j 



(43) 



The channels which preserve the maximally mixed state are called 
bistochastic channels. The structure of one-qubit bistochastic channels is 
discussed in more detail in Section [4. II 



1.8 Correlation matrices 

A general measurement process is described in quantum mechanics by 
operators forming a positive operator valued measure (POVM). Products 
of matrices K % ' K % representing the POVM are positive and determine 
the identity resolution, 11 = K l . During the measurement of 

a quantum state p the output p,: = ^^f^pf occurs with probabilities 

Pi = Tr K % pK 1<[ . The identity resolution guarantees that Y^ii=iP 1 = 1- 
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The outcomes of a quantum measurement are not perfectly distin- 
guishable, unless different POVM operators project on orthogonal sub- 
spaces, IC^ICK^K' = SijK^K*. Probability distribution of the out- 
come states does not contain any information on indistinguishability of 
outcomes. Therefore, a better characterization of the measurement pro- 
cess is given by the following correlation matrix a with entries 

o-ij =Tr RtpK*, i.j 1 fc. (44) 

Its diagonal contains the probabilities of measurement outputs, while the 
off-diagonal entries are related to probabilities that the state % has been 
determined by the measurement as the state j. The correlation matrix 
depends on both, the measured state and the measurement process. 

The operators K l , satisfying X/»=i K^K* = 1, can also be treated 
as Kraus operators (|30p characterizing the quantum channel, $(p) = 
5^i=i K % pK % ^. In such an interpretation of operators K % , the correla- 
tion matrix (I44|) is equivalent to the state of environment given by the 
output of the complementary channel $(p) specified in Eq. (|36|l . 

The entropy of the state a produced by a complementary channel $ 
is called the exchange entropy, since, if the initial states of the system 
and the environment are pure, then S(a) is equal to the entropy which is 
gained by both the state and the environment [29] . If the initial state is 
maximally mixed, p = p» = irl, where N is the dimensionality of p, the 
entropy of the output of the complementary channel is equal to the map 
entropy S rmap ($) [46] (see also discussion in Section [2. 1.1 J |. 

S map (*) = -^D*log(lW), (45) 

where the dynamical matrix D® is given by Eq. (|25[) . This entropy is 
equal to zero if $ represents any unitary transformation. It attains the 
largest value log 2N for completely depolarizing channel which transform 
any state into the maximally mixed state. Therefore the map entropy can 
characterize the decoherence caused by the channel. 

Due to the polar decomposition of an arbitrary non normal operator 
X = HU, we can write K'p 1 ^ 2 = hiUi, where hi is a Hermitian matrix 
and Ui is unitary. One can observe that h% = K'pK^ — Pipt.. Therefore 
the entries of the correlation matrix (|44[) can be written as: 

a l3 = Tr K'pK^ = pj pj Tr pf UdJ]p) . (46) 

As noticed above, the correlation matrix characterizing the quantum 
measurement can be equivalently treated as the state of an environment 
after evolution given by a quantum channel. The following section indi- 
cates a third possible interpretation of the correlation matrix a. It can be 
formally treated as a Gram matrix of purifications of mixed states pi. 

Purification of a given state pi £ Mn is given by a pure state 
(see Eq. (21)), 

Tn |*,) = pi. (47) 
The purification [SPj) of given state pi can be written explicitly, 

N 

l*i> =^2(Ui® yfpi) |r>® |<#.), (48) 

r=l 
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where {|$.)}r^=i are eigenvectors of pi. Notice that a purification of a 
given state pi is not unique. The degree of freedom is introduced by the 
unitary transformation Ui. Moreover, any purification of given state pi 
can be given by such a form. Since eigenvectors of pi denoted by 1 4>V) 
form an orthonormal basis in the Hilbert space, a unitary transformation 
Vi can transform it into the canonical basis {[r)}^. The purification (148p 
can be described as 

N 

|*i> = { Ui ® V^i) \r) ® \r} ■ (49) 

r=l 

The overlap between two purifications of states pi and pj emerging from 
a POVM measurement is given by 

| 2 = | (m\ (U]Ui ® V/vW/^i) |m) | 2 , (50) 

where \m) = \r) ® jr). For any two operators A and _B the following 
relation holds, (m| A®B \m) = Tr B 62J . Hence the overlap (|50[) reads 

K* J |* i )| a = |'nw A /«>/£1 a , (5i) 

where the unitary matrix W = VillJ UjVj . Therefore the matrix elements 
of a (|46p are equal to the scalar product of purifications of respective 
mixed states pi and pj as follows a%j = ^/PiPj^ j^i) ■ 

1.8.1 Gram matrices and correlation matrices 

In previous chapter it was shown that the correlation matrix can by de- 
fined by the set of purifications of states emerging from the quantum 
measurement. Therefore, the correlation matrix can be identified with 
the normalized Gram matrix of the purifications. 

The Gram matrix is an useful tool in many fields. It can receive a 
geometrical interpretation, as it consists of the overlaps of normalized 
vectors. If vectors are real the determinant of their Gram matrix defines 
the volume of the parallelogram spanned by the vectors |63II64| . The 
Gram matrix of the evolving pure state is analyzed in |65j . The spectrum 
of this matrix can determine whether the evolution is regular or chaotic. 

The Gram matrix a, 

°a = y/PiPj(^i\ipj) (52) 
has the same eigenvalues as 

P = £>hMM- (53) 

i 

The proof of this fact [66] uses the pure state, 

\4>}=^VFiWi}®\ei), (54) 
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where states \ei) form the set of orthogonal vectors. Since the state (|54[) 
is pure, its complementary partial traces equal to (|52|) and (|53|) have the 
same entropy 

s ([Vm^(^j)h) = s fe>|iMW*M ■ (55) 

The entropy of the Gram matrix (|52[) can be used in quantum information 
theory to describe the ability of compression of quantum information |67j . 
The authors of [67] describe the fact that it is possible to enlarge the 
information transmitted by means of set of states which are pairwise less 
orthogonal and thus more indistinguishable. This fact encourages us to 
consider global properties of quantum ensemble which, sometimes, are not 
reduced to joint effects of each pair considered separately. In Chapter [3] 
some efforts will be made to define the quantity characterizing fidelity 
between three states. 



1.9 Kraus operators constructed for an ensemble 
of states 

The previous section concerns the ensembles £ — {pi,Pi}i = i formed by 
the outputs of a given quantum channel and a given input state. In 
the following section it will be shown that for any ensemble £ the suitable 
Kraus operators K % can be constructed and the corresponding initial state 
p can be found. 

Initial state is constructed from the states of the ensemble by taking 

k 

P=52piU} Pi Ui, (56) 

i=l 

where the unitary matrices Ui are arbitrary. The Kraus operators con- 
structed for ensemble £ and unitaries Ui are defined by 

K^^plUn (57) 

Vp 

Notice that K 1 pK^ — pipi and the Hermitian conjugation, K'^ = -j=U} ^JpTpi 
Due to the choice of p in (156[) the identity resolution holds, 

^2 K H K i = ^2 p .J_ u j p . u .J_ = h (58) 

In the special case of k = 2 states in an ensemble, by choosing 

U 2 = Ui 1 ^_ ^pTV&, (59) 

one obtains a\2 equal to square root fidelity between states p\ and p2, as 
follows y / F(p 1 ,p 2 ) = Tr \J y/pip2y/pi- 

In consequence of the above considerations one can say that the en- 
semble emerging from POVM measurement can be arbitrary and for any 
ensemble £ we can construct the set of operators K l and the corresponding 
initial state p. 
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1.10 Quantum fidelity 

An important problem in the theory of probability is how to distin- 
guish between two probability distributions. The so called fidelity is a 
quantity used for this purpose. Assume that P = (pi,P2, ...,pjv) and 
Q = (qii, q2, gjv) are two probability distributions. The fidelity between 
p and q is defined as, 

F(P,Q) = (j2VP^J ■ (60) 

This function has several properties: 

• it is real, 

• positive, F(P,Q) > 0, 

• symmetric, F(P,Q) = F(Q,P), 

• smaller or equal to unity, F(P, Q) < 1. 

• equal to one if and only if two distributions are the same, 
(F(P,Q) = 1)^(P = Q). 

These properties are shared by fidelities defined for quantum states given 
below. 

Quantum counterpart of the fidelity for the pure states \4>i) £ Hn and 
\4>2) G Hn is given by the overlap 

F(|0i>,|<M) = |^i|02)| 2 . (61) 

A probability distribution can be considered as a diagonal density matrix. 
Generalization of two formulas (|60[1 and (|61[) for arbitrary mixed states 
Pi G M. n and p 2 G Mn is given by 



F(jn,P2) = (Tr^PWpi") ■ (62) 

To show a relation to previous definitions of fidelity consider two com- 
muting quantum states. They can be given, in the same basis, as pi = 
2~Zf ri I*) an d Pi = 2- Si N) ('I- Hence the fidelity between them reads 



(Tr y/Vplp^/piy 



= Tr, 



N 



This gives a relation between fidelity between mixed quantum states (|62[) 
and fidelity of probability distributions which are composed by the eigen- 
values of the states (|60|) . Consider now pure states, |$i) , ^2) G Hn^Hn 
such that the partial trace over the first subspace reads, Tri — pi. 

There exists a relation between formula (|62[) for fidelity between two mixed 
states and overlaps of their purifications. 

Theorem 2 (Uhlmann [62]). Consider two quantum states pi and p2 and 
their purifications |*P '1) and I'lPa). Then 



(TrJVfrPzVP^) =max|(*i | * 2 )| 2 , (64) 
where the maximization is taken over all purifications |^i) of the state pi. 
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Proof. The proof starts from purification formula (|49[) . 

i**} = (Ui ® V^^i) N . (65) 

where |m) is an unnormalized vector, \m) = YliLi \ r ) ® 10- The overlap 
of two purifications (|50|l is given by 

\{y 3 \<K t )\ 2 = \TrW^/pJ^-\ 2 , (66) 

where the unitary matrix W = ViUjlIjVj . The maximization over pu- 
rifications is equivalent to maximization over the unitary matrix W. An 
inequality | Tr AB\ < \\A\\ Tr \B\ provides the required lower bound 

| Tr WJp-^pl\ 2 < (Tr | Jp]Jpl\ f . (67) 

The upper bound is attained by the unitary matrix W' equal to the 
unitary part of the polar decomposition of ^J~pjyfpl- This finishes the 
proof. □ 

1.10.1 Geometrical interpretation of fidelity 

Consider two one-qubit states in the Bloch representation 1 140( 1 . 



p x = i (id (68) 

Pv = i (id (69) 

where S is the vector of Pauli matrices (1391) . Fidelity of the pair of states 
p x and p y reads 



F(p*,Pv) = ~{i + £-y + y/i z PFvi - llyll 2 )- (70) 

If the states p x and p y are both pure then ||£|| = ||y|| — 1 and the fidelity 
can be given by 

F{px,Py) = COS 2 -, (71) 

where the angle a is formed by two Bloch vectors which represent the 
pure states p x and p y at the Bloch sphere. One can use this statement 
to define the angle between two states as a function of the fidelity. The 
generalization of such an angle for arbitrary two mixed states is given by 



A(p!,p 2 ) ■■= arccos y / F(p 1 ,p 2 ). (72) 

It was proved |68| that such an angle satisfies the axioms of a distance 
and leads to a metric. 



1.11 Mutual information 

The goal of quantum information is to efficiently apply quantum resources 
for information processing. Consider the following situation. A sender 
transmits the letters of the message from the set X — {oi, 02, — , The 
letters occur with probabilities pi, where i = 1, k. The message is trans- 
mitted by a communication channel, which can be noisy and can change 
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some of the letters. The receiver performs a measurement and obtains 
outputs Y with a possibly different probability distribution. According to 
the Shannon information theory Q] the amount of information contained 
in the message characterized by probability distribution pi is given by the 
entropy H(X) = — ^2-p% log Pi- Entropy describes the average amount of 
digits per letter necessary to transmit the message characterized by this 
probability distribution in an optimal encoding scheme. 

The receiver knowing the letters Y has only a part of information 
contained in the original message X. The information which Y and X 
have in common is characterized by the mutual information H(X : Y) 
defined by 

H(X : Y) = H{X) + H(Y) - H(X, Y), (73) 

where H(X, Y) is the Shannon entropy of the joint probability distribution 
of the pairs of letters, one from X and one from Y. 

The errors caused by a channel can be perfectly corrected if the mutual 
information is equal to the entropy of the initial probability distribution. 
Otherwise the mutual information is bounded by the entropy of an initial 
distribution [8], 

H(X : Y) < H(X). (74) 
Following properties of the mutual information hold [8J: 

• Mutual information does not change H(X : Y, Z) = H(X : Y) if the 
system Z is uncorrelated with Y. 

• Mutual information does not increase if any process is made on each 
part, H[X : Y) > H(X' : Y'), where prime denotes the states after 
the transformation. 

• If part of a system is discarded the mutual information decreases 
H(X :Y,Z)> H(X : Z). 

Mutual information can also be defined for quantum composite sys- 
tems in terms of the von Neumann entropy . The definition is analogous 
to (73J): 

S(p P : p Q ) = S(p P ) + S(pq) - S(ppq), (75) 

where states of subsystems are given by partial traces, for example, pp = 
TtqPpq. Mutual information S(pp : pq) for quantum states satisfies 
properties analogous to these listed above for the classical mutual infor- 
mation H(X,Y). 

1.12 Holevo quantity 

Holevo x quantity (Holevo information) of the ensemble £ = {g^pij^Lj 
is defined by the formula 

/ k \ k 

X({^pi}) = S ^ q lPl - J2 ( 76 ) 
\i=l / i=l 

It plays an important role in quantum information theory. As the bound 
on the mutual information |7J, Holevo quantity is related to fundamental 
restriction on the information achievable from measurement allowed by 
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quantum mechanics. It directly reflexes these features of quantum me- 
chanics which distinguishes this theory from classical physics. In classical 
information theory the mutual information between the sender and the 
receiver is bounded only by the Shannon entropy of the probability dis- 
tribution describing the original message. In the case of an ideal channel 
between two parts the mutual information is equal to the upper bound. 
In quantum case, even without any noise present during the transmission 
process, the mutual information is restricted by the Holevo quantity which 
is smaller than the entropy associated with the original message, unless 
the states used to encode the message are orthogonal. 

The theorem of Holevo [7] is presented below together with its proof. 

Theorem 3 (Holevo). Let {pi}*L x be a set of quantum states produced 
with probabilities pi from the distribution P. Outcomes of a POVM mea- 
surement performed on these states are encoded into symbols with probabil- 
ities qj from probability distribution Q. Whichever measurement is done, 
the accessible mutual information is bounded from above, 

/ k \ k 

\i=l / i = l 

Proof. Consider a three partite state, where its parts are denoted by the 
letters P, Q and M 

UpQM = ^ Pi \i) P i ®\0}{0\. (78) 

i 

Three parts of the system P, Q and M can be associated with the prepara- 
tion state, quantum systems, and the measurement apparatus respectively. 
The state ljpqm describes the quantum system before the measurement, 
since the state of the apparatus is independent on the quantum states. 

Assume that the state lopqm is subjected to the quantum operation 
acting on the subsystem QM as follows, <E>(p ® |0) (Oj) = £). K j pK^ ® 
\j) The Kraus operators of this quantum operation form a POVM 
measurement since K j = \. The state after this measurement is 

given by 

up-q-m' = Y,P> I*) (*l ® K 3 p t K 3 ^ ® \j) (j\ . (79) 

y 

Properties of the mutual information listed in section fl . 1 1 1 imply the 
key inequality of the proof: 

S{ljp : ojq) > S{lu p , : uM'). (80) 

To prove inequality (|77[) it is enough to calculate the quantities occurring 
in (|80[) for the state (|78[) and (|T9[) respectively. Since lopq = Tr mijJpqm = 
X^i Pi I*) (*l ® P» j the left hand side of (|80[) is given by 

k 

S(u;p : ojq) = S(uj p ) + S(u Q ) - S(u PQ ) = S(p') - ^^(p,), (81) 

1=1 

where p' — ^tPiPi- This is the Holevo quantity which does not depend 
on the measurement operators K % . To compute the right hand side of 
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(|80[) . S(ujpi : cjm'), consider a state (|79[) . The observation that p(x,y) — 
p x p(y\x)=p x TrK yt K y p x leads to 



S(u P ,:u M ,)=H(P:Q), (82) 

where Q = {q v } v and q y = Tr K v p' K v * . This is the mutual information 
between the probability distributions describing the outcomes of the mea- 
surement and the original message. That finishes the proof of the Holevo 
bound on the mutual information of message encoded into quantum sys- 
tems. □ 

Above theorem is one of the most important applications of the Holevo 
quantity. Quantum information theory uses also the Holevo quantity % 
to define channel capacity. There exist several definitions of quantum 
capacity of a channel depending on whether the entanglement between 
the input states is allowed or not. In the case that quantum states in a 
message are not entangled the Holevo capacity of channel $ is defined by 



Ch(&) — max 

£={p»,«}tU 



/ k \ k 

\i=l / 8=1 



(83) 



The Holevo quantity x(£)> which can be interpreted as the Holevo capacity 
of the identity channel, bounds the capacity Ch for any channel |5j: 

Ch < x{£)- (84) 

Yet another application of the Holevo quantity concerns the ensembles 
of quantum states. Formula (|76|) can be given by the average relative 
entropy 

k / k \ / k \ k 

^Pi D I Pi^PiPi ) = S I ^Pipi ) - ^PiSipi), (85) 

i = l \ j = l / = l / i = l 

where the relative entropy is defined as D(pi,p2) = Tr pi (log pi — logp2). 
It defines an average divergence of every state from the average state. 
Average (|85|) is known as the quantum Jensen Shannon divergence [69 . 
Its classical version, for probability measures, is considered in 70 . From 
mathematical point of view, the Holevo quantity can be treated as a quan- 
tity which characterizes the concavity of the entropy function. 

The Holevo information will be the main object considered in Part [H] 
of this thesis. 
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Part II 

Bounds on the Holevo 
quantity 

2 Holevo quantity and the correlation ma- 
trix 

In the following chapters several inequalities for the Holevo information 
(Holevo quantity) will be given. It is well-known [5] that the Shannon 
entropy of the probability vector P — {pi, pk} is an upper bound for 
the Holevo quantity of an ensemble £ — {pt, /Oj}f =1 : 

X {£) < H(P). 

Since the Holevo quantity forms a bound on accessible mutual information, 
the difference between entropy of probability vector H(P) and the Holevo 
quantity specifies how the chosen set of density matrices differs from the 
ideal code, which can be decoded perfectly by the receiver. The upper 
bound on the Holevo quantity can be used for estimating this difference. 
One of the estimation for the Holevo quantity is presented in the following 
section. 

As discussed in Section lTTSl the correlation matrix a can be equivalently 
interpreted in several ways. If the set of the Kraus operators K 1 defines a 
quantum channel, $(p) = 53i=i K 1, pK*\ the correlation matrix a charac- 
terizes the output state of the complementary channel, a = $(p), or the 
state of the environment after the quantum operation. As mentioned in 
Section ri.8.11 a defines also the Gram matrix of purifications of the states 
{pi}i =1 . The entropy S(a) is related to the exchange entropy or the en- 
tropy which the environment gains during a quantum operation provided 
the initial state of the environment is pure. In the following analysis a 
quantum channel $(p) = J3 i K' pK l * is treated as a device preparing an 
ensemble of quantum states £ = {pi, pi}i-i, where 

Pi = Triryrt, and Pi = J£j^L , (86) 

The described situation is illustrated in Fig. [1] 

Independently of the interpretation of the Kraus operators K r the 
following theorem proved in [49 holds. 

Theorem 4. Let ^2 i=l K^ K % — 1 be the identity decomposition and 
p an arbitrary quantum state. Define the probability distribution pi — 
Tx K x pK 1 ^ and a set of density matrices pi — -jf^^-it • The Holevo 
quantity x({P»>Pi}) * s bounded by the entropy of the correlation matrix, 

^ = El=iTr#y^t N)(j | ; 

k k 
1=1 1=1 
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Figure 1: A quantum channel $ represents a device preparing the ensemble 
of quantum states £ = {pi,Pi}f =1 - The average of this ensemble is denoted as 
p' = &(p) — J2i=i K l pK % ^ . The complementary channel <£> transforms an initial 
state p into the state a of the environment. 



where H(P) is the Shannon entropy of the probability distribution P = 
{pi,—,Pk}- 

Proof. The right hand side of the inequality: S(a) < H(P), is a con- 
sequence of the majorization theorem, see e.g. [15] . Since the probabil- 
ity vector P forms a diagonal of a correlation matrix, we have S(a) < 
5(diag((r)) = H(P). The left hand side of the inequality (|87p is proved 
due to the strong subadditivity of the von Neumann entropy [60] . The 
multipartite state CJ123 is constructed in such a way that entropies of its 
partial traces are related to specific terms of 11871) . 

The multipartite state W123 is constructed by using an isometry F\<f>) = 
J2i=i H) ® I*) ® The state CJ123 = FpF f is given explicitly by the 

formula 

fc 

^123 = FpF^ = J2 \i){j\ ® ® K l P K ] \ (88) 

States of the subsystems cji are given by partial traces over the remaining 
subsystems, for example, ui = Tr23Wi23 and so on. 

Let us introduce the following notation Aij = K' pK^ . In this no- 
tation the quantities from the Theorem U take the form pi = Tr An and 
pi = An I Pi. Notice that 

S(wia) = Sip), (89) 
fc 

Siu 3 ) = S(Y t PiPi)- (90) 

i=l 

Moreover 

fc fc fc 

-^ Pl S(pO = ^ Tr An log A<j - ^ Tr(A«) log Tr(A«) 

i— 1 i — 1 i— 1 

= - SV23). (91) 

The strong subadditivity relation in the form which is used most fre- 
quently 

SV23) + Sim) < S(lj 12 ) + 5(w 23 ) (92) 
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does not lead to the desired form (|87[) . However, due to the purification 
procedure and the fact that a partial trace of a pure state has the same en- 
tropy as the complementary partial trace, inequality (|92l) can be rewritten 
in an alternative form |21| : 

S(ws) + S(ui) < S(uia) + S(was). (93) 

This inequality applied to the partial traces of the state (|88|) proves 
Theorem H □ 

For an ensemble of pure states pi = \tpi){ipi\, the left hand side of 
(|87[) consists of the term 5 , (^ i pi|i/'i}(V'!|) only. The correlation matrix 
o in the case of pure states is given by the Gram matrix. Due to the 
simple observation (155 fl , the left inequality (I87p is saturated in case of any 
ensemble £ consisting of pure states only. 

Using a different method an inequality analogous to Theorem [4] has 
been recently proved in [71] for the case of infinite dimension. It can 
be also found in |72| in context of quantum cryptography. The authors 
analyse there the security of a cryptographic key created by using so called 
'private qubits'. In such a setup an inequality analogous to (|87[) appears 
as a bound on the information of the eavesdropper. 



2.1 Other inequalities for the Holevo quantity 

Methods similar to that used to prove Theorem [4] can be applied to prove 
other useful bounds. 

Proposition 1. Consider a POVM measurement characterized by opera- 
tors X/i=i K^K 1 = II which define the outcome states, pi — ^^i^ K i\ an d 
their probabilities, pt = Tr K l pK*'. The average entropy of the output 
states is smaller than entropy of the initial state, 

k 

^>S(pO < S{p). (94) 

8=1 



Proof. Due to the fact that the transformation F in Eq. (|88[) is an isom- 
etry, the three-partite state CJ123 has the same nonzero spectrum as the 
initial state p. Hence CJ123 and p have the same entropy. Due to equality 
(|91[) and the Araki-Lieb inequality 76J: 

5(wi) - S(w 23 ) < SV123), (95) 
one completes the proof of Proposition [T] □ 

Note that concavity of entropy implies also another inequality 53i=i Pi&{p 
S(p') = SQ2i=iPiPi)- Proposition [T] has been known before |77| as the 
quantum information gain. 

Definition of the channel capacity (|83|) encourages one to consider 
bounds on the Holevo quantity for the concatenation of two quantum 
operations. Treating the probabilities pi and states pi as the outputs 
from the first channel one can replace maximization over £ = {pi,pi}^ 1 
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in (|83[) by maximization over the initial state p and the quantum oper- 
ation $1. The strategy similar to that used in Theorem U allows us to 
prove the following relations. 

Proposition 2. Consider two quantum operations: $i(p) = Ej=i KlP^l 
and $ 2 (» = E£i KipKp. Define Pl = Tr K\pK^ and p % = ■ 
The following inequality holds: 

ki ki 

S{$ 2 o $!(/,)) - ^p,5($ 2 ( Pl )) < S($i(p)) - X>S(/x)- (96) 

i=l i=l 



Proof. Let us consider the four-partite state: 

fci k 2 

^1234 =EE K> 01 ® ® ^^P^iff, (97) 

where |nn) = |n) ® |n), and the strong subadditivity relation in the form 
SV; 24 ) + S(t4) < S(w 14 ) + (98) 

Notice that 

S(^) = S($ 2 o$i(p)), 
- S(uj' 24 ) = -Y.iPiS{®2{Pi)), 
S(u' 14 ) = S(E- 2 3=1 K> 0'] ® Ktei(j>)K*) = S($i(p)). 

The third equality is due to the fact that an isometry, F2\<j>) = E;=i I 1 ') ® 
^2 10) j does not change the nonzero part of spectrum. This property is 
also used to justify the following equation 

fci 



Substituting these quantities to the strong subadditivity relation (|98|> we 
finish the proof. □ 

Inequality [96] is known [8] as the property that the Holevo quantity 
decreases under a quantum operation x(PmP0 > x(P») ^K/ 3 *))- 

Consider notation used in the proof of Proposition [2] Concavity of the 
entropy gives 

k\ k\ k 2 fci fc2 

X>s($ 2 (pi)) =5> iS (S>^) - EX>® S (^)- ( 10 °) 

z — 1 i — 1 3—1 i = 1 j = 1 



where = t^k^k^k^k^ and P robabilities = Tri^^iA^i^- 
Using Theorem 2] and concavity of entropy (|100|) one proves: 
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Proposition 3. Consider two quantum operations: $i(p) = Y2i=i K\pK^ 
and $ 2 (p) = E^i KipKp- Define Pi = Tr K\pK? and p; = ^rj^ ■ 
The following inequality holds: 

fci 

S($ 2 o $i(p)) - ^piS^aCpi)) < S(an), (101) 

where the output of the complementary channel to $2 £3> $1 is denoted as 
a H = $ 2 o $i(p). 

2.1.1 Some consequences 

This section provides three applications of theorems proved in Sections [2] 
and l2.ll One of them concerns the coherent information. This quantity is 
defined for a given quantum operation $ and an initial state p as follows 
[73] 

/ coh (*,p) = 5(*(p))-5(*(p)), (102) 

where 4>(p) is the output state of the channel complementary to $. To 
some extent, coherent information in quantum information theory plays 
a similar role to mutual information in classical information theory. It is 
known ^8 that I C oh(&, p) < S(p). That is a relation similar to (|74[) . More- 
over, it has been shown that only if I co h{&,p) = S(p) the process $ can 
be perfectly reversed. In this case the perfect quantum error correction is 
possible [73] . The coherent information is also used to define the quantum 
capacity of a quantum channel |74| 

C*q($) = max J coh ($,p). (103) 

p 



The definition of the coherent information (|102[) can be formulated 
alternatively [73] by means of an extended quantum operation $ ® id 
acting on a purification \tf>) 6 W2 <8> H3 of an initial state, p = Tr3[i/>}(^i|. 
This fact is justified as follows. The purification of p determines as well 
the purification Q123 £ Hi ® Hi £S> Hz of the state uj £ Hi <8> H2 in (|29|) . 



fim = C/12 ®id 3 (|l)(l|i ® |V) (^123)^2 ®id 3 - (104) 
The partial trace over the environment (subspace Hi) reads 

«2 3 = [$®id](h/>)(V>|). (105) 

It has the same entropy as the partial trace over the second and third 
subspace, Q,i = a, which is a state of environment after evolution, 

S(a) = 5([*®id](|V)(^|)), (106) 

and S(o-) = 5($(p)). 

Coherent information (|102[) can be written as 

/«*(*, p) = S(Tr 3 «23) - 5(n 23 ). (107) 
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The classical counterpart of the coherent information can be defined by 
using the Shannon entropy instead of the von Neumann entropy and prob- 
ability vectors instead of density matrices in Eq. (|107p . The classical 
coherent information is always negative, since the entropy of a joint prob- 
ability distribution cannot be smaller than its marginal distribution. 

Inequalities proved in Theorem U and Proposition [T] together provide 
the following bound on the coherent information, 

k 

W*,p)<X>s(*)<s(p), ( 108 ) 

1=1 

where pi = Tr K'pK^ and pi — K 1 pK'^ /p, are defined by Kraus repre- 
sentations of the channel, <&(p) = X^=i K'pK'^ . The equality between 
coherent information and the entropy of initial state S(p) guarantees that 
<3> is reversible. Inequality (|108[) implies a similar, weaker statement: only 
if the following equality holds Xli=i Pi^(Pi) = S(p), the quantum opera- 
tion <E> can be reversed. 

Another consequence of inequalities proved in Section T2 . 1 1 concerns the 
so called degradable channels. These channels are considered in quantum 
information theory in the context of their capacity [42] . A channel $de 9 
is called degradable if there exists a channel ^ such that ^ o &deg = <&de 9 . 
Substituting the degradable channel <I?i = <&de 9 and the additional channel 
$2 = ^ to inequality in Proposition [2] one obtains a lower bound for the 
average entropy of 4'(pi), where pi are output states from the channel 

®deg, 

k k 

< Y.P^) - W***,p) < X>5(¥(/>0), (109) 

i=l i=l 

where I co h{<&, p) = S^deg(p)) — S(&deg(p)) ■ The left inequality is due to 
inequality (|108p . Therefore Proposition [2] provides some characterization 
of the channel which is associated with a degradable channel. 

The third application of propositions from Section [2. II is given as fol- 
lows. The Jamiolkowski isomorphism |34| gives a representation of a quan- 
tum map $ which acts on N dimensional system by a density matrix on 
the extended space of size iV 2 . This state can be written as: 

a* = [id®$](j^ + )(<£ + |), (110) 

where \<t> + ) = 4j>^j=i I*) ® I*) * s ^ ne maximally entangled state. A 
rescaled state D$ = Na<s> is called the dynamical matrix. In the special 
case, if the initial state is maximally mixed, p = 4=-l, the entropy of 
the correlation matrix a written in (|106|) is equal to the entropy of the 
dynamical matrix. 

A quantum map $ can by defined using its Kraus representation (|30[) . 
Since the Kraus representation is not unique |15| . one can associate many 
different correlation matrices with a given quantum operation <£> depend- 
ing on both, the initial state and the set of Kraus operators. However the 
entropy of the dynamical matrix D<s, is invariant under different decom- 
positions. This entropy characterizes the quantum operation and is called 
the entropy of a map [49], denoted by S map ($) as defined in Eq. (|45|) . 
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Due to Theorem [4] the entropy of a map has the following interpre- 
tation. It determines an upper bound on the Holevo quantity (|76|) for a 
POVM measurement defined by the Kraus operators of <E> if the initial 
state is maximally mixed p = p* = -^1. Moreover, the entropy of a map 
is an upper bound for the Holevo quantity for POVM given by any set of 
Kraus operators {K 1 ' K 1 } which realize the same quantum operation $, 

max X ({ Pi =Tr K%K*, Pl = **£f^ }) < S(*), (HI) 

where p' = $(p) = J2i=i K* pK^ . 

Proposition [3] provides also an alternative lower bound for the entropy 
of composition of two quantum maps given by Theorem 3 in |46j . The 
inequality for the entropy of composition of two maps can be now stated 
as 

k 

< MAX{S($ 2 o $i (/>*)) - ^2 Pi S(^ 2 (pi)), S($i) + A | < 5($ 2 ° *i), 

i=l 

(112) 

where A = S($ 2 o $i(p,)) - S($i(p.)) and $(p) = Y, i= iPiPi- The 
lower bound proved in our earlier paper |46| could be smaller than 0. The 
improved bound is always greater than due to concavity of entropy. 

2.2 Discussion on the Lindblad inequality 

Lindblad [75] proved an inequality which relates the von Neumann en- 
tropy of a state p, its image p' — <E>(p) = Yl-=iPiPi an( l the entropy of 
the correlation matrix a equal to the output state of the complementary 
channel a = $(p), 

\S(p) - S(p)\ < S(a) < S(p') + S(p). (113) 

Another two Lindblad inequalities are obtained by permuting the states 
p, p and a in this formula. The proof of Lindblad proceeds in a similar 
way to the proof of Theorem 0] It involves a bi-partite auxiliary state 
w" = J2i =i N)01 ® KipKj, where the identity S(p) = S(lj") is due to 
an isometry similar to F in (|88[) . The Araki-Lieb inequality |76| . \S(pi) — 
S(p2)\ < S(pi2) applied to u>" proves the left hand side inequality of 
(|113|) . while the subadditivity relation S(pi 2 ) < S(pi) + S(p2) applied to 
lj" proves the right hand side inequality of (|113[) . 
Inequalities from Theorem [3] and Proposition Q] 

k 

S{p) - ^PiSiPi) < S(<t), (114) 

k 

^PiS(pi) < S(p) (115) 

i=l 

use a three-partite auxiliary state u = Y^i j=i N*)0il ®K%pK^. As in the 
case of the Lindblad inequality (|113[) . the identity S(p) = S(lj) holds due 
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to isometry. The strong subadditivity relation applied to oj proves inequal- 
ity (|114p . while the Araki-Lieb inequality applied for uj proves inequality 
(|115|) . Notice that an extension of the auxiliary state and application of 
the strong subadditivity relation allows one to use the average entropy to 
new inequalities for interesting quantities: the entropy of the initial state, 
the entropy of the output state of a quantum channel p' = $(p) and the 
entropy of the output state of the complementary channel ( 3?(p). 

In the case S(p') > S(p) (e.g. for any bistochastic operations) the 
result (|114p gives a better lower constraints for S(a) than the Lindblad 
bound (I113[) . In this case 

fc 

S(p') - S(p) < S(p') - ^PiSfa) < S(a), (116) 

i=i 

due to Prop. [T] However, if S(p') < S(p) the result of Lindblad can 
be more precise depending on the values of S(p), S(p') and the average 
entropy Yli=iPi^(Pi)- in consequence, due to Lindblad inequality (| 1 13f) 
and the inequality (|111|) one obtains another lower bound for the entropy 
of a map: 

MAx{log(JV)-Sr(*(p.)), max X (p, = TrK'p*K'\ Pi = T f^ p f^ t ) } ^ ^ map ($), 

"(117) 

where p> = $(p) = £*_i K 'p K ^ = Y*=xViQi- 

2.3 Inequalities for other entropies 

Inequality (|87[) uses the strong subadditivity relation in the form (|93f) 
which is a specific feature of the von Neumann entropy. Relation (|93[) can 
be equivalently formulated in terms of relative von Neumann entropies. 
The relative von Neumann entropy D(pi, p%) is defined as follows 

D{pi, P 2) = Trpi[logpi - lng(pa)] (H8) 

and is finite for p2 G supp(pi), otherwise it becomes infinite. 

Monotonicity of relative entropy states that for any three-partite quan- 
tum state U1123 and its partial traces the following inequality holds: 

D(l)23,OJ2 ® CJs) < D(bJi23, Wl2 ® W3). (H9) 

It is an important and nontrivial property of the von Neumann entropy 
[60] . |78| . Monotonicity of the von Neumann entropy (|119|> rewritten using 
the definition (|118|) leads to the strong subadditivity relation: 

SV23) + S{u) Z ) < SV13) + S(u 23 ). (120) 

Complementary partial traces of any multipartite pure state have the 
same entropy. This fact can be applied to purifications of W123. Therefore, 
relation (|120[) is equivalent to (|93[) which can be applied to the specific 
three-partite state (|88|) 

fc 

^123 = J2 K>01 ® 1001 ® K iP K ] (121) 
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and used to prove the upper bound on the Holevo quantity in terms of a 
correlation matrix \ < Hence, inequality (|87p is a consequence of 

the monotonicity of the relative von Neumann entropy. 

Monotonicity of entropy holds also for some generalized entropies e.g. 
Tsallis entropies of order < a < 1 ,79J or Renyi entropies of order 
< q < 2 [80J. Direct generalization of \ < •S'C' 7 ) I s n01, so easv i since 
the key step in the proof was the strong subadditivity form (|93[) . In 
case of generalized entropies such a form cannot be obtained from the 
monotonicity of relative entropy. 

The Holevo quantity can be expressed by the relative entropy. Con- 
sider the state (|121|) and the notation: K 1 pK^ = Pipi, and ^2^ =1 Pipi = 
p' . The relative entropy reads: 

D(u>23,UJ2 ®w 3 ) = (122) 

= Tr UJ23 log 0J23 — Tr ui2s log tJ2 — Tr UJ23 log W3 (123) 
fe fc 



= ^Tr ^ pi log pi pi -^pilogpi -Trp'logp' (124) 

i=l i=l 
fe 

= ^ Pi Tr pi log pi - Tr p log p (125) 

i=l 

fe fe 

= S(pO-£p*S(*)=£p^(«,p')=X- (126) 

i=l i=l 

The equality between the Holevo quantity and relative entropy holds also 
for the Tsallis entropies of any order q 

T a (p) = I ^-[l-Tr P °'], (127) 

where the relative Tsallis entropy of order a is defined as [79| 

^(p 1)P2 ) = -l--[l-Trp^- a ]. (128) 

It is now possible to compute the Tsallis-like generalized relative entropy 
between a bipartite state U23 and the product of its partial traces 
which leads to the generalized Holevo quantity %a- H one considers the 
state (fT2T|) 

DZ{u2-i,u2 ®w 3 ) = — l — -[l-Trw? 3 (wa®W3) 1_a, l (129) 

Ot — 1 L J 

= ^[i-E^fepO^v 1 ^] ( 13 °) 

i=l 

= $>^-r(l -Tr pfp' 1 "") (131) 
z — ' a — 1 

i=l 
fe 

= J2^ D -^'p"^ = ^- ( 132 ) 
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In a similar way we can work with the Renyi entropy Sf(p) = log[Tr p a ]. 
The corresponding relative Renyi entropy reads [81] 

D*(pi,P2) = -^jlogTrlplp 1 2 - q ] (133) 

and the Renyi-Holevo quantity is given by 

xf = -^jlogTr(J2p^) 1/q - (134) 

i 

Equality between the generalized Renyi-Holevo quantity (|134[) and the 
Renyi relative entropy (|133p holds if relative entropy concerns partial 
traces of (|121|) and the state p" = (Yl< i PiPi) 1 ^ q as follows 

X f = D?(u> 23 ,U2®p"). (135) 

The Holevo quantity (|135|) is smaller than Df (^23,^2 0J3) |81| . 

The monotonicity of relative entropy for three considered types of 
generalized entropies: von Neumann entropy, Tsallis entropy of order 
< a < 1 and Renyi entropy of order < q < 2 gives 

X< D(Lu 123 ,a(g>p'), (136) 
< Dl{ui2z,(J<& p), (137) 
xf <D?(u 123 ,<T®p). (138) 

These relations state that the Holevo quantity is bounded by the relative 
entropy between the joint state of the quantum system and its environ- 
ment and the states of these subsystems taken separately. 

In case of von Neumann entropy, inequality (|136|) can be written ex- 
plicitly as 

X<S(a) + S(p')-S(p). (139) 

Notice that p is an initial state and S(p) = 5(^123) due to isometry 
transformation, F : p — !> 1^123- Relation (|139|) joints entropies of the 
initial state, the final state, the state of the environment and the Holevo 
quantity in a single formula. Inequality (|139[) which can be rewritten as 

k 

S(p) < S(a) + Y,PiS(Pi) (WO) 

i = l 

gives a finer bound than that provided by the Lindblad inequality: S(p) < 
S(a) + S(p'). Inequality ()139p can be written as x < S(a) + Y , where 
y| = \S(p') — S(p)\ < S(a), due to one of the Lindblad inequalities. In 
some cases this inequality confines the relation (1871) . 

2.4 Searching for the optimal bound 

The state a can be defined for a triple consisting of a probability distri- 
bution, set of k density matrices of size iV and a set of k unitary matrices, 
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{pi, pi, Ui}i =1 . Every triple (p,, pi, Ui) defines uniquely the pure state \ipi) 
which is the purification of state pi as follows 

JV 

1^) =Y^ (IJ *®^P lV ^\ e r)® M ( 141 ) 
r=l 

as shown in (149[) . The Holevo quantity depends only on £ = {pi,Pi}i=i. 
Therefore, Theorem |4] can be reformulated as follows: 

Theorem 5. For any ensemble {pi, pt,Ui}i =1 the Holevo quantity is 
bounded by the entropy of the correlation matrix a minimized over all 
unitary matrices Ui 

x({PuPi}) = S(^2piPi)-^2piS(pi) < minS'(cr) = mm Pi\ipi) (V>i|), 

i—1 i= 1 ' i— 1 

(142) 

where \tpi) = Ylr=AUi<8v r PiVi)\e r )®\er) andoij = y/piPj Tr ^fPi^fpjU]Ui . 

The last equality of (|142l) holds since the correlation matrix a can be 
represented as the Gram matrix of purifications of pi. It is known that 
for any Gram matrix equality (|55[) holds. 

Finding minimization of S(a) over unitaries is not an easy problem 
in general. In the following chapter the problem will be solved for the 
ensemble of k — 2 states, and the solution is written in terms of square 
root of the fidelity between both states. A conjecture that the matrix of 
the square roots of fidelities also bounds the Holevo quantity for ensembles 
of k = 3 states will be formulated and some weaker bounds will be proved 
in the next section. 



2.4.1 Optimal bound for two matrices 

The tightest upper bound on the Holevo quantity occurring in Theorem 
[5]is obtained by taking minimum of S(a) over the set of unitaries. This is 
equivalent to the POVM which minimizes the correlation matrix among 
all POVM which give the same output states. For two output states pi 
and p2 occurring with probabilities (A, 1 — A) the correlation matrix is 
given by 

VA(1 - A) Tr y/fiy/p&lUi 1 - A 

(143) 

Its entropy is the lowest, if the absolute values of the off-diagonal ele- 
ments are the largest. As has been shown in Eq. (|67[) the expression 
Tr ^/~p\yfp2UlU\ attains its maximum over unitary matrices at the value 

(144) 

where for brevity we use F12 instead of F(pi,p2). This quantity is equal 
to the square root fidelity (|62p . Therefore the correlation matrix of the 
smallest entropy can be rewritten in terms of the square root fidelity, 



A_ y/X{l - A)v / 7\? 

\A(1- A)v / 7TI 1- A 



(145) 
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2.5 Jensen Shannon Divergence 



Minimal entropy of the correlation matrix characterizing an ensemble of 
two density matrices is related to the distance between them in the set 
of density matrices. If the probability distribution in (|145[) is uniform, 
A = 1/2, the square root of the von Neumann entropy of a m in forms a 
metric |53j . It is called the entropic distance De(pi, P2) 



DE(pl,p2) = \J S(amin), 



0~ m i n — Z 



1 V f (pi>p*) 



(146) 

Inequality ()142|) provides the relation between this metric and another 
one defined by means of the Jensen-Shannon Divergence. The Jensen- 
Shannon Divergence JSD({a v P„}) has been initially defined [69], [82] 
as the divergence of classical probability distributions P v occurring with 
probabilities a v 

JSD{{a„P v }) = H ( ctuPu) - a » H ( p ») = a » H ( p » \\P) ( 147 ) 



where H(P) denotes the Shannon entropy of the probability distribution 
P, H(P V \\P) is the relative entropy between P v and P, while the average 
probability distribution reads P = a v P v . 



De(PQ) 




D E (P,Q)-D T (RQ) 




Figure 2: a) The entropic distance De (1146)) for two 2-point probability distri- 
butions P — (p,l—p) and Q = (q, 1 — q). b) The difference between the entropic 
distance De and the transmission distance Dt (| 1*491) - 

The square root of the Jensen-Shannon divergence between two prob- 
ability distributions P and Q, 

JSD(P\\Q) = \H{P\\M) + \H{P\\M), (148) 

where M = ^(P + Q), forms a metric in the set of classical probability 
distributions [82] . |83| called the transmission distance Dt(P,Q), 

D t (P, Q) = y/JSD(P\\Q). (149) 

A probability distribution can be considered as a diagonal density matrix. 
Therefore, Eq. (|142[) in Theorem [5] demonstrates a relation between func- 
tions of two distances in the set of diagonal density matrices. Fig. [5] and 
Fig. [3] shows the comparison between these two distances for exemplary 
probability distributions. 
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A quantum counterpart of the Jensen-Shannon divergence, in fact co- 
inciding with the Holevo quantity, was also considered |69| . |82) . Inequal- 
ity (|142[> provides thus an upper bound on the quantum Jensen-Shannon 
divergence. 




Figure 3: a) The entropic distance De (|146ll and the transmission distance Dt 
(|149p for two probability distributions, P = (i, i , ^) and Q — (<7i, 92, 93) which 
is arbitrary distribution of dimension 3 represented by a point in the simplex - 
the base of the figure, b) The difference between the entropic distance De and 
the transition distance Dt for the same distributions P and Q. 



3 Conjecture on three—fidelity matrix 



The minimization problem for the entropy of the correlation matrix (|143|) 
has been solved for an ensemble consisting of k = 2 quantum states. In 
this case the solution is given by the square root fidelity matrix. In the 
case of k = 3 states in the ensemble the optimization over the set of 
three unitary matrices is more difficult. Our numerical tests support the 
following conjecture, which is a generalization of the bound found for the 
case of k = 2. 

Conjecture 1. For an ensemble of k = 3 quantum states, {pi,pi}^ 1 the 
entropy of the square root fidelity matrix Gij = ^/PiPj \/F(pi, pj) gives the 
upper bound on the Holevo quantity, 



x({Pi,Pi}) < s 



Pi tJp~VP2V F12 

VP2PlVf21 P2 
^PSPlVFil ^P3P2\/F 3 2 



VP1P3V F13 

^/pip-iVF^ 

P3 



(150) 



where fidelity between two quantum states reads Fij = F(pt, pj) = (Tr y/^/PiPj \fPi) ' 

It has been shown |84| . |52| that the matrix G containing square root 
fidelities is positively semi-defined for k = 3. However, the square root 
fidelity matrix is in general not positive for k > 3. Numerical tests provide 
several counterexamples for positivity of G for k > 3, even in case of an 
ensemble of pure states. Note that the matrix G is not a special case of 
the correlation matrix a, which is positive by construction. 
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Theorem [4] implies that Conjecture \T\ holds for ensembles containing 
three pure states. Inequality (|87[) is in this case saturated as discussed in 
section[2] Square root fidelity matrix G is obtained from the Gram matrix 
of given pure states by taking modulus of its matrix entries. Taking 
modulus of entries of a positive 3x3 matrix does not change neither 
the trace nor the determinant of the matrix. Only the second symmetric 
polynomial of the eigenvalues is growing. Since the entropy is a monotonic 
increasing function of the second symmetric polynomial |67| . the entropy 
of the square root fidelity matrix G is larger than the entropy of the Gram 
matrix and therefore it is also larger than the Holevo quantity. 



3.1 A strategy of searching for a proof of the con- 
jecture 

The proof of Theorem [3] consist of two steps. In the first step one has to 
find suitable multipartite state. In the second step the strong subadditiv- 
ity relation of entropy has to be applied for the constructed multipartite 
state. The same strategy will be used searching for the proof of Conjecture 
[T]or for proving other weaker inequalities. 

For the purpose of obtaining the Holevo quantity from suitable terms 
of the strong subadditivity relation, the multipartite state uj should have 
a few features: 

• it is a block matrix which is positive, 

• blocks on the diagonal should contain states pi multiplied by prob- 
abilities pi, 

• traces of off-diagonal blocks should give square root fidelities, or 
some smaller numbers if one aims to obtain a weaker bound. 

The following matrix satisfies above conditions, 



X 



Pi Pi 





o 
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o 








* 








o 
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P2p2 
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* 
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o 








P3P3. 



(151) 



where in place of * one can put any matrix, provided the matrix X remains 
positive. If in place of * one substitutes zeros, the strong subadditivity 
relation implies the known formula that x({Pi> Pi}) < S({Pi})- Examples 
presented in the next section use described strategy to prove some entropic 
inequalities for the Holevo quantity. 

The main problem is to find a suitable positive block matrix. In order 
to check positivity the Schur complement method [85] is very useful. 

Lemma 1 (Schur). Assume that A is invertible and positive matrix, then 



X 



A 



B 

C 



(152) 
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is positive if and only ifS — C — A 1 B is positive semi-definite: 

A>0 => (X>0 <=> S > 0). (153) 
The matrix S is called the Schur complement. 



3.1.1 Three density matrices of an arbitrary dimension 

The strategy mentioned in the previous section will be used to prove the 
following 

Proposition 4. For a three states ensemble {pi, Pi}i=i,2,3 the following 
bound for the Holevo quantity \ holds 



x(Pi,Pi) < s 



where b > 2. 



Pi ^PiP2y/F 12 /b ^pip-WFrs/b 

^/P2PlVT2l/b P2 y^P2pl^F^/b 

_y/P3P2\/F3l/b yJP3P2VF~32~/b f>3 



(154) 



Proof. It will be assumed that considered density matrices {pi}^ 1 are in- 
vertible. After 106 the square root of the product of two density matrices 
yffxr will be defined as follows: 



pa : 



(155) 



In this notation the fidelity between two states pi and pj can be written 
as: 

/ i \ 2 

\2 



F(pi,Pj) 



- r I 1/2 1/2 

Tryp/ Pjp/ 



(TrVpIpJ) 



(156) 



Formula (|156[) can be generalized for non-invertible matrices |52| . 

One can use the Schur complement Lemma[T]to prove positivity of the 
block matrix: 

pi ^fpTpi 
VP2P1 P2 



X 



(157) 



In this case the matrices A and S, which enter the Lemma [TJ take the 
form: A = pi, assume that it is invertible, and S = p2 — \JPzPiP\ ypiP2. 
Notice that 

P2 - S = % /p 2 piprV/ 9 iP2Pipr 1 (158) 
= V 'P2P1V V^PiPi -1 = P2, (159) 
S = and X > 0. Hence the 



therefore in the case of matrix (11571 
following matrix Y is also positive: 

i 



Y 



\pi 




|yp2P 











\pl 









2\/P3Pl 





\/pTp2 



\p2 


















IV^pTpl 









\P2 








\yfPlP~3~ 







|P3 










\jp2p~i 


\pz 
(160) 
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Using strong subadditivity as described in section 13.11 to the multipar- 
tite state TT2Y extended by some rows and columns of zeros, one proves 
inequality (|154|) for b — 2. To prove relation l)154p for b > 2 a small 
modification of matrix (|160[) is needed. The off-diagonal elements can be 
multiplied by the number < r < 1 without changing the positivity of 
the block matrix. □ 



3.1.2 Three density matrices of dimension 2 

Proposition [4] can be amended for the case of 2 x 2 by decreasing the 
parameter b to the value at least v3- 

Proposition 5. For an ensemble of three states of size two, {pi,pi}f = i 
one has 



x(Pi,Pi) < s 



Pi ^/PiP2VFi2/b yJpipW Fr.i/b 

VP2P1V-P21/6 P2 v / P2P3v r F23/b 

_^p-iP2\?F [n /b ^JpzP2^fFv2lb p 3 



with b > v3. 

Proof. The main task in the proof is to show that the block matrix 



(161) 



Y = 



pxpi v / PiP2v / piP2/b y/pipl-^/Pipl/b 

^P2ply^P2pl/b P2p2 ^p2Ply^p2P3/b 

_y / plP2yfpTpl/b ^P3P2y^P3p2/b p 3 p3 



(162) 



is positive for b > Vs as well as the analogous matrix enlarged by adding 
rows and columns of zeros in order to have a matrix of the form (|15ip . The 
Schur complement method described in section [3~T1 will be used, where: 



A = 



1 

p\p\ 



C : 



P2P2 ^/P2P3y / P'2p:i/b 

VP2P3V 'p-iP-ijb p 3 p 3 



(163) 



B = 





^/plpl^Plpl/b ^pTply^pTpl/b 



B = 



^/PlP2^/P2Pl/b 

^JPiP3^P3pi/b 

(164) 

Due to the fact that A is positive one needs to prove the positivity of 
S = C-B^A~ 1 B: 



P2P2(1 - p-; 

/V2P3{yJP3P2/b - y/P3PlPi 



IP\P2 



/P2P3(v 'P2p 3 /b - s]p2p\Px X \JP\P3lb 2 ) 
P3P3(1 - pr) 

(165) 



To prove positivity of (|162[) the Schur complement S should be positive. 
One can apply the Schur complement Lemma second time to the matrix 
S. Positivity condition required by Lemma [T] enforces that 



b 2 (b -3)p!+by>0, 



(166) 



where y = ^/pTp2p 2 1 \/P2p 3 p s 1 y / P3~pT + h.c. For 2x2 matrices one can 
assume without lost of generality that -^=y-j= > 0. It is so because the 

matrix \fp~i~ 1 ^/ pip2P2 X \/ p2P3p 3 1 p 3 p\ \fp~\~ X is a unitary matrix and 
its determinant is equal to 1, therefore its eigenvalues are two conjugate 
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numbers. The matrix -^=y— ==, which consists of sum of the unitary 
matrix and its conjugation, is proportional to identity. If it is negative 
one can change ^Jpxpi into — y ' p\pi and % /p 2 Pi into — ^/p 2 pi in (|162p . 
Transformation changing the sign does not act on the final result because 
off-diagonal blocks do not take part in forming the Holevo quantity and 
in the case of 3 x 3 matrices we can take modulus of each element of the 
matrix without changing its positivity. 

Let us take y — in the positivity condition (|166p . This condition 
implies b > \/3. Knowing that (|I62[) is a positive matrix, the rest of the 
proof of (|I6I[) goes like in section [3TT1 □ 



3.1.3 Fidelity matrix for one qubit states 



In previous section some bounds on the Holevo quantity were established. 
These bounds are weaker than the bound postulated by Conjecture [TJ 
since decreasing the off-diagonal elements of a matrix one increases its 
entropy. In previous proposition the square root fidelities were divided 
by numbers greater than 1. In the following section the squares of the 
off-diagonal elements of the matrix G in (|150|) will be taken. For such 
modified matrices the following proposition holds for an arbitrary number 
of k states in the ensemble. 

Proposition 6. Consider the ensemble {pi,pi}^ 1 of arbitrary number 
k of one-qubit states and their probabilities. The Holevo information 
xdPiiPi}) * s bounded by the entropy of the auxiliary state ? which acts in 
the k - dimensional Hilbert space, 

x({pt,Pi})<S{<i), (167) 
where ^ = ^/pppj (Tr yjplp]) 2 = y/PiPj F(pi,pj). 

Proof. A positive block matrix W is constructed in the following way: 



W = 



Mi 
M K 





• a4" 


. 





. 






(168) 



where Mi = y/pi(Ai, B%) are block vectors of size 2x4 and Ai — pi and 
Bi = y^detpTll are sub-blocks of size 2x2. The blocks of the block matrix 
W read 

Wij = y/plpj{p,pj + ^/detpipji). (169) 
This formula can be compared with an expression for the square root of 
any 2x2 positive matrix X 



TrVX 

Therefore the block matrix (I168[) is given by 

Wij = ^PiPj^/PzpjTi-y/PiPj- 



(170) 



(171) 



The matrix W is positive by construction. Partial trace of this matrix 
gives matrix of fidelities (without square root). The rest of the proof of 
Proposition [5] goes in analogy to proofs analysed in Section [XT] □ 
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This proposition holds for one-qubit states only since we applied rela- 
tion (|170|) . which holds for matrices of dimension d = 2. 

The fidelity matrix ?y = \JPiPj (Tr -*JPiPj) 2 = ^/PiVj F*ij is n °t positive 
for a general k and general dimensionality of pi. However the fidelity ma- 
trix is positive and bounds the Holevo quantity in the case of an ensemble 
containing an arbitrary number of pure quantum states of an arbitrary 
dimension. This is shown in the following proposition. 

Proposition 7. Let {\tfj)} be a set of vectors, then 

x({Pi,Pi})<S(F), (172) 
where = */plp]\{<Pi\<Pj)\ 2 - 

Proof. Introduce a complex conjugation tp Tp by taking complex con- 
jugations of all coordinates of the state in a given basis. Hence for any 
choice of tp, tp one has 

= <V#>. (173) 
The matrix F2 := [F(pi, pj) 2 ]ij can be rewritten as 

[\(fi\fj}\ 2 h = [(V»IVi>(Vj|^i>]«J 

= [{<Pi\<Pi)tjPiW])]ii (174) 

= [((<Pi\ ® (<Pi\)(\<Pi) ® 

This last matrix is the Gram matrix of the set of product states {\tpj) ® 
|^7)}j =1 and therefore is positively defined. 

The next part of the proof continues according to the scheme presented 
in Section r3.ll We use the multipartite state 

w = VPiPi (Vi I fi) I**) J I ® IVtXVjl- (175) 
y 

Its positivity is shown by taking the partial trace of the Gram matrix 

^ = ^2 l«)0'il ® I* 3 *) ® l^iXwl ® (&/|- (176) 

a 

The proof is completed by considering partial traces of the state oj and 
using the strong subadditivity relation. □ 



3.1.4 Special case of the correlation matrix 

The previous propositions use the strategy from the proof of Theorem 
2] and apply it to positive block matrices which are not necessary re- 
lated to the correlation matrices. Construction of multipartite states al- 
lows one to obtain the matrices containing fidelities after a partial trace. 
The following section deals again with the correlation matrices aij = 
y/PiPj Tr y^y/pjt/j Ui . Since the Holevo quantity does not depend on 
unitaries Ui, these matrices can be chosen in such a way that the three- 
diagonal of a consists of the square fidelity matrices, Oij = \jF{fii, pj), 
where \i — j\ < 1. This construction is used in the following proposition. 
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Proposition 8. Consider an ensemble {pj,Pi}i = i consisting of arbitrary 
number k of invertible states of an arbitrary dimension. The Holevo 
information x({Pi> Pi}) * s bounded by the exchange entropy S(a), 

X {ipi,Pi})<S(a), (177) 

where the correlation matrix a is given by: 

on = p 4 , (178) 

o%i = ■s/piPjO^-s/plpl), iff |i-j|=l, (179) 

and the upper off-diagonal matrix elements, where (j — i) > 1, read: 

1 1 1 



o-ij = x 7~lT, Tr VPjPj -i 



\/P3-lPj-2 



y/p^Fi, (180) 



pj-l pj-2 pi+1 

while lower off diagonal satisfy = oji . 

The matrix a has a layered structure presented here for k = 4, 



Pi 








0" 




" 


/l2 





" 




" 





f (2) f (3)"| 
J 13 J14 





P2 
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/ai 





/23 





+ 








o /£> 








P3 










/32 





/34 


f m 

J31 

















P4_ 










/« 


_ 




A3) 
J 41 


f(2) 
J 42 



(181) 



with entries of this matrix equal to cry specified in Proposition [8] 

Proof. Consider a correlation matrix: 

o-ij = ■s/plplTr. s /p'i^/pJUjUi (182) 

where unitaries C/i are chosen in such a way that elements at i±i are square 
root fidelities: ^jFT7±i = Tr ^/pT7h±i- Hence 



TV — V' IV 



where Vj_ 1 a is the unitary matrix from the polar decomposition, 



s/~Pi\[Pi = \VP~iVP~J\ v i-J = V p] /2 Pip] /2 Vi,j. 
Here the Hermitian conjugated unitary matrix V t reads: 

»••» 7777 . /7T V ^ 



1/2 1/2 
Pi Pi ■ 



(183) 



(184) 



(185) 



The first unitary Ui can be chosen arbitrarily. The recurrence relation 
(11831) allows one to obtain formula (11801). □ 



To analyse properties of the matrix a consider, for example, the matrix 
element o~\z- 

0"i3 = \fpYP~3 Tr Jpi^fpi V\U\ 

= y/plpl Tr VpT^pI Vis Ul Ui 
= VPiP3 Tr ^fpl^fpi Vl , V^a . 



(186) 
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Using Eq. fl 185 p one obtains 



/ rp , , J- J- / 1/2 1/2 1 1 / 1/2 1 

c-13 = VPiPs Tr Jpl^pi—— —— y p 2 ' p 3 p 2 —^=—^=V Pi P2p! 



that gives the matrix element CT13 of (|18QJ) - The assumption that the ma- 
trices are invertible is used in (118511 where the unitary matrix of the polar 
decomposition of ^/fH^fpJ is given explicitly. However, the same strat- 
egy of the proof leads to analogous proposition involving non-invertible 
matrices. Only the equations ()185|) and (|180|) are changed in this case. 

3.1.5 Hierarchy of estimations 

One can compare average values of entropies from Conjecture Q] and 
Propositions [4] [6] and [8] The average values are situated on the scale 
in which the Holevo quantity is set to and the entropy S(P) of proba- 
bility distribution is set to unity. The variable s *fy_ x is used, where x is 
replaced by the entropy of respective state. The standard deviations are 
also computed. The probability distributions are generated according to 
the Dirichlet measure, while the set of k = 3 density matrices is chosen 
randomly according to the Hilbert-Schmidt measure |86| on the set of 
density matrices of size 2. 

• <X>=0 

• < Sfid >= 0.176 ± 0.065, where Sfid corresponds to the entropy 
from Conjecture [1] 

• < Siayered >= 0.193 ± 0.087, where Siayered corresponds to the 
entropy from Proposition [8] for k = 3 states in the ensemble. 

• < Sfid 2 >= 0.37 ± 0.13, where Sf id 2 corresponds to the entropy 
from Proposition [5] for k — 3 states in the ensemble. 

• < Sfid/b >= 0.750 ±0.015, where Sfid/t corresponds to the entropy 
from Proposition U 

• < S(P) >= 1. 

For an ensemble of k — 3 one-qubit states Conjecture[T]is the strongest, 
as it gives on average the lowest bound, while among the statements 
proved in Propositions 2] E] and [5] the tightest bound (on average) is pro- 
vided by Proposition [8] 

3.2 Fidelity bound on the Holevo quantity for a 
special class of states 

Although, Conjecture [TJ has been confirmed in several numerical tests, it 
has been proved so far for the set of pure states (Section [3]) only. The aim 
of the following section is to prove that the square root fidelity matrix 




(187) 
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bounds the Holevo quantity for a restricted set of states. It will be shown 
that for one-qubit states among which two are pure and one is mixed and 
for the uniform probability distribution, {|, |, i}, Conjecture \T\ holds. 

Proposition 9. Consider k — 3 one-qubit states pi among which two 
are pure pi = \4>i)(<f>i\, P2 = \4>2)(4>2\, and the state p3 is mixed. The 
square root fidelity matrix G for these states and the uniform distribution 
P = {|, |, |} bounds the Holevo quantity, 



x(fo,ft}) < s ~ 




(188) 



(Try7p7p~7pT) 2 . 

The proof goes as follows. First proper parameters characterizing three 
states will be chosen. After that the formulas for the left and right side of 
inequality 1)1880 . which are functions of two variables only, will be given. 
The fact that one of these functions is greater than the other is shown 
graphically. 

Notice that the left hand side of Eq. (|188)l depends only on the lengths 
of the Bloch vectors which represent the mixed state p$ and the average 
state p = ^(pi + p2 + p'i) inside the Bloch ball. The same average p can 
be realized by many triples {pi,P2,P3} where pi,P2 are pure and p3 is 
mixed of given length of the Bloch vector. The family of such triples is 
parametrized by two numbers a and (5 as shown in Fig. [4] The points 
B, D, E denote the following states: B — > pz which is mixed, D — > p 1 = 
|<^i}(0i| and E — > p2 = \4>2){4>2\, while A — > p represents the average state. 
The vector OA of length a denotes the Bloch vector of the average state p, 
the vector OB of length 6 characterizes the mixed state P3. The position 
of the vector OB with respect to OA can be parametrized by an angle a. 
These two vectors, OB and OA, determine, but not uniquely, two pure 
states from the same triple characterized by OD and OE. Equivalently 
one can rotate the vectors OD and OE by an angle /3 around the axis 
OC and obtain pure states denoted by F and G. The ratio \AB\ : \AC\ is 
equal to 2 : 1 because in this case the average A is the barycenter of three 
points B, D and E or a triple B, F and G. The method of obtaining 
the points C, D, E, F and G, when a, b and a are given, is presented in 
Appendix 1. Given a pair of parameters (a, b) distinguishes the family of 
triples {|0i)(</>i|, \4>2){<f>2\, pi} characterized by two angles a and /3. The 
range of a is given by condition \OC\ < 1, it is 



iV9a 2 -6fecos(a)a + fe 2 < 1 

< Q < TV, 



while the range of j3 is (0, 7r). Left hand side of Eq. ()188p depends only 
on the lengths a and b, and is independent of the concrete realization of 
the triple. Therefore to prove (|188[1 for given a and b one has to find 
minimum of the entropy of the square root fidelity matrix over all triples 
parametrized by the angles a and f3. 
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Figure 4: The Bloch representation of the three states {pi — \<fii)(4>i\, P2 — 
1^2) (<?^2 1, P3} and the parameters used in the proof of Proposition[!)]are presented 
schematically in the Bloch ball. Two angles (a, (3) characterize all possible triples 
{pi, p2, /03} if parameters (a, b) are fixed. 



The entropy of the square-root fidelity matrix defined Gij = ^/PiPj \jFi 
in Eq. (|188|l is a function of roots of the characteristic polynomial: 

(i-A) 3 +p(i-A)+q = 0, (190) 

where 

p = -(F12 + F13 +F 2S )/9 (191) 

q = 2VF12F13F23/27. (192) 

The parameter p determines the second symmetric polynomial S2 of eigen- 
values of the square root fidelity matrix G 

82 = \Yp-~ F *i)- ( 193 ) 

i<j 

The roots of equation (|190l) are equal to: 

X k = i + 2^cos [(I arccos (|^X) + ,f )], (194) 

where fe = 1, 3. 

The entropy of the square root fidelity matrix is a function of p and 
q, which determine the second symmetric polynomial of eigenvalues (|193[) 
and the third symmetric polynomial is in this case equal to the determi- 
nant of the 3x3 matrix Gij = y/PiPj \J F\j ■ The von Neumann entropy is 
a monotonically increasing function of all symmetric polynomials of eigen- 
values [67] . The parameter q is a function of (o, b, a, ft), while parameter 
p depends only on a and b which is shown in following lemma: 

Lemma 2. For any triple of two pure and one mixed state of an arbitrary 
dimension the sum of fidelities depends only on the purity of the mixed 
state and the barycenter of the ensemble. 
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Proof. Denote by p the barycenter of a mixed state p and two pure states, 

I0i>, \<h), 

p = lp+l\<h){4>i\ + l\4») {<l*\- ( 195 ) 

The purity of p is given by 



Trp 



i ( Tr p 2 + 2 + 2 ((f)! | p {fa) + 2 (0 2 | p]02> + 2 |(0i| 2 )| 2 ) . (196) 



After reordering the terms one gets 



where 



Frz + F 13 + F 23 = i(9Trp 2 - Trp 2 - 2), 



(197) 



Fl2 = | (01 1 02)| , 
F23 — (02 1 P |02) , 
-Fl3 = (01 1 P|0l) • 



(198) 

(199) 
(200) 



Since Tr p 2 — | ( 1 + a 2 ) and Tr p 2 = | ( 1 + 6 2 ) , the parameter p defined in 
(1191[) does not depend on the angles a and /3. This completes the proof 
of Lemma [5] □ 

The parameter p and the second symmetric polynomial (|193p does not 
depend on the angles a and /3. Therefore, for given a and b, the entropy 
of the square root fidelity matrix attains its minimum over a and j3 for 
minimal value of the determinant of G, since the entropy is an increasing 
function of the determinant. The determinant is given by 



det 




_L_ 

27 



1 + 2V-F12-F13F23 - (-F12 + F13 + F 23 ) 
(201) 



It is the smallest for the smallest value of the parameter q which is the 
function 1)1921) of the off-diagonal elements of the matrix. During compu- 
tations of the minimal value of q another lemma will be useful: 

Lemma 3. Among triples of one-qubit states which realize the same 
barycenter, where one state is mixed of a given purity and two others 
are pure, the product of three pairwise fidelities is the smallest if three 
states and the average lie on the plane containing the great circle of the 
Bloch ball, i.e. ft — 0. 

Proof. The function f(a,b,a,/3) — F12F13F31 is given explicitly in Ap- 
pendix 2 based on Appendix 1. For given a, b and a this function has 
minimum only at j3 = and equivalently for /3 = it. □ 

In consequence, searching for the minimum of the entropy of the square 
root fidelity matrix we can restrict our attention to the case /3 = 0. In 
fact, for our purpose it suffices to take the specific value of a which is 
shown in the following lemma. 
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Lemma 4. Among triples of one-qubit states which realize the same 
barycenter, in which one state is mixed of given purity and two others 
are pure, the product of three pairwise fidelities is the smallest when two 
pure states are symmetric with respect to the mixed state i.e. /3 = and 
a = or a = ir. 

Proof. The function fo(a, b, a, /3 — 0) = F13F13F31 is given directly in 
Appendix 1. It has only one minimum at a — but in certain cases, 
depending on a and b, the value on the edge of variable range, i.e. at 
a = or a = n is smaller. □ 

3.2.1 Proof of the fidelity bound 

To prove inequality (|188p the smallest entropy of the square root fidelity 
matrix for three states consistent with the left hand side of this inequality 
should be found. Entropy is a function of four parameters, (a,b,a, 0). 
The left hand side of (|188|) . which is the Holevo quantity depends on two 
parameters (a, b) as follows 



For given parameters a and b lemmas 1, 2 and 3 allows us to find specific 
a and /3 for which minimization of right hand side of (|I88[) is obtained. 
One can fix a = or a — ir and /3 — 0. That means, that minimal entropy 
of the square root fidelity G over the angles is obtained if the three states 
{pi = \(j>i}(4>i\, p2 = \4>2)(<f>2\, P3} are lying on the great circle and the two 
pure states are symmetric with respect to the mixed state. In this case 
the matrix G is characterized by two parameters, F = F12 = F23 and b. 
Here F is the fidelity between the pure state pi and the mixed state P3 
whereas b characterize the length of the Bloch vector of the mixed state 
P3. The matrix G reads 



where F is a function of b, such that F(b) = |(1 — fee), and c is the 
length of the Bloch vector representing the barycenter of two pure states 
pi and p2- The fidelity F is equal to 1/2 if b tends to 0. The parameter 
c determines also the projection of the Bloch vector of the pure state pi 
on the Bloch vector of the mixed state P3. The absolute value \c\ is equal 
to the square root fidelity between the two pure states. The range of 
variables are < b < 1 and |(1 - b) < F < + b). 

Considered case is shown in Fig. [5] There are two surfaces - func- 
tions of two parameters F and b. The lower surface represents the Holevo 
quantity \, and the upper surface denotes the entropy of the square root 
fidelity matrix (|203[) . The surface S(G) lies always above x an d is com- 
posed of two smooth functions characterizing cases in which all vectors 
lay on the same semicircle or pure states and the mixed state belong to 
the opposite semicircles. 




(202) 




(203) 
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Fig. [5] suggests that in the case of three pure states, 6=1, laying on 
the same semicircle the inequality is saturated, \ — S(G). In this case, 
F > 1/2, the rank of the square root fidelity matrix is equal to 2, and the 
nonzero eigenvalues are (l±a)/2, where a = (4F— 1)/3 is the length of the 
Bloch vector of the average state p. In general we have a = |(6 + 2 2F b ~ 1 )■ 
In case of x = S(G) the Holevo quantity is equal to the entropy of the 
average state p. This finishes the proof of Proposition [9] 




Figure 5: Evidence in favour of Proposition |9] The Holevo quantity as function 
of two variables: fidelity F between the pure state pi and the mixed state p3, and 
the length b of the Bloch vector characterizing the state pz- The upper surface 
representing the square root fidelity matrix G is composed of two smooths parts. 
Every circle represents schematically the Bloch ball with exemplary positions 
of Bloch vectors characterizing three states {pi — pi = \4>2)(4>2\, P3} of 

the ensemble. 
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Part III 

Minimal output entropy and 
map entropy 

4 Entropies for one-qubit channels 



The question on additivity of the channel capacity is one of the most 
interesting problems in quantum information theory |40| . Shor showed 
[39] that this problem has several equivalent formulations. One of them 
concerns the minimal output entropy, 



In the case of one-qubit channel the minimal output entropy is the entropy 
of a state characterized by point on the ellipsoid, which is the image of 
the Bloch sphere, the closest to this sphere. The pure state which is 
transformed into a state of the minimal entropy is called minimizer. 

For any setup in which minimal output entropy is additive the quan- 
tum channel capacity is additive as well. Additivity implies that an en- 
tangled state cannot increase capacity of two channels with respect to 
the sum of their capacities taken separately. The additivity conjecture 
can also be formulated as a statement that capacity of two channels is 
minimized for a product state. 

The conjecture was confirmed in many special cases. For instance, 
additivity holds, if one of the channels is arbitrary and the second one is: 
bistochastic one-qubit map [87] . a unitary transformation [40| . generalized 
depolarizing channel [41] . entanglement breaking channel 88J, very noisy 
map [89| and others. A useful review on this subject was written by 
Holevo [90| . Different strategies of proving the additivity conjecture are 
analyzed there. For a recent relation on the additivity conjecture see 
also [15] . 

Also counterexamples to the additivity conjecture have been found. 
One of them was presented by Hastings [16) . He found the lower bound 
for the output entropy of some channels when the input was a product 
state. Next he estimated the output entropy for a maximally entangled 
input. Due to such estimations it was shown that the entangled state 
decreases channel capacity below the value achievable for product states. 

The proof of Hastings used pairs of complementary channels. His 
argument was not constructive and works in high dimensional spaces. 
Counterexamples for the additivity hypothesis are also studied in |17| . 

It is still an open question, whether the additivity holds for an arbitrary 
one-qubit channel. Originally, the hypothesis on additivity of minimal 
output entropy was formulated for the von Neumann entropy. One of 
the approaches to the problem uses a one-parameter family of entropies, 
called Renyi entropies characterized by a parameter q, 
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minS($(»). 
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Calculations are sometimes easier when the Renyi entropies are consid- 
ered. The quantity S q tends to the von Neumann entropy in the limit 
q — > 1. Additivity of the minimal output Renyi entropy has been proved 
only in some range of the parameter q depending on the channels consid- 
ered namnHzi- 

Although the Renyi entropy is sometimes computationally more fea- 
sible, finding minimum over entire set of quantum states is still a hard 
problem. One of the ideas how to omit this difficulty tries to use some 
relations between minimal output entropy and other quantities which are 
easier to calculate. In the following chapter the Renyi entropy of a map 
(the map entropy) will be used to estimate the minimal output entropy. 
Map entropy (entropy of a map) is defined by the entropy of the Choi- 
Jamiolkowski state (|28[) corresponding to the map. This quantity is easy 
to obtain. Numerical tests presented in Fig. [JJ [TO] show that there is 
no simple functional relation between the map entropy and the minimal 
output entropy. Nevertheless being aware of the structure of the set of 
quantum maps projected on the plane (S^ la ' p , S™ ln ) can be useful. Knowl- 
edge of entropies of maps at the boundaries of the allowed set can be used 
to estimate the minimal output entropy by the entropy of the map. 

4.1 Structure of the set of Pauli channels 

Quantum channels which preserve the maximally mixed state are called 
bistochastic. All bistochastic one-qubit channels can be represented as a 
convex combination of the identity matrix cro = II and three Pauli matrices 
CTi=l,2,3 (ESJ) 

3 3 

®P(P) = ^PiViPVi, = !) ^iPi > °- ( 206 ) 

i=0 i=0 

Bistochastic one-qubit quantum operations are thus called Pauli channels. 
The structure of the set of all Pauli channels forms a regular tetrahedron 
A3 as shown in Fig. |BJi. There are many channels characterized by the 
points of tetrahedron which can be obtained from other channels follow- 
ing a unitary transformation. Our considerations are often restricted to 
the asymmetric tetrahedron K (see Fig. 0>) which is a subset of A3. 
All maps in A3 can be obtained from channels of K by concatenation 
these channels with unitary transformations. The set K is formed by the 
convex combination of four vectors p from H206|) , A = (0,0,0,0), B = 
(1/2,1/2,0,0), C = (1/3,1/3,1/3,0), and D = (1/4,1/4,1/4,1/4). 
Extremal lines of the asymmetric tetrahedron correspond to the following 
families of maps: AB - dephasing channels, BD - classical bistochas- 
tic maps, AD and CD - depolarizing channels. The families mentioned 
above are also shown in Fig. [7J which presents boundaries of the set of all 
one-qubit bistochastic channels projected onto the plane (<S'™ ap , S™ In ). A 
following proposition proved in |51| characterizes this projection. 

Proposition 10. Extremal lines of asymmetric tetrahedron correspond 
to boundaries of the set of all bistochastic one-qubit maps on the plot 

(sr p ,sr n )- 
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(0,0,0,1) 




Figure 6: The structure of one-qubit bistochastic quantum operations corre- 
sponds to the regular tetrahedron A 3 . This figure is spanned by four extremal 
vectors p from formula (|206p . Symmetries of the tetrahedron allow us to dis- 
tinguish the asymmetric set K inside A3 . Any vector p characterizing a Pauli 
channel can be obtained by permutation of elements of vectors from K . 




Figure 7: Lines AB, BD and AD, which correspond to the edges of asymmetric 
tetrahedron K form the boundaries of the entire set of Pauli matrices projected 
on the plane (S , 2 nap ,5™ n ). 
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4.2 Depolarizing channels 

Fig. [9] and Fig. 1 101 show the projection of the Pauli channels on the plane 
(<S'™ ap , 5 1 ™ 111 ) with parameter q different than 2. Comparison of these 
figures with Fig. [7] shows that the structure of the set of channels on the 
plane (S'™ ap , S™ ln ) is the simplest in case of the Renyi entropy of order 
q — 2. In this case, the depolarizing channels form one of the edges of the 
set of all quantum one-qubit maps projected onto the plane. Indeed the 
following theorem proved in |51| confirms the special role of depolarizing 
channels in the set of all quantum channels acting on states of arbitrary 
dimension N. 

Theorem 6. Depolarizing channels have the smallest map Renyi entropy 
S'™ ap among all channels with the same minimal output Renyi 



The family of depolarizing channels is represented in the plane (S™ ap , S 1 ™ 111 " 1 
by the continuous line on the entire range of S™ ap . The minimal output 
entropy of a depolarizing channel Ajv acting on M m is the following func- 
tion of the map entropy 

S 2 mm (sr P (Aiv)) = - log ( N+1 ) • (207) 

This is a monotonously increasing function from to logiV. Therefore 
the following theorem holds. 

Theorem 7. Depolarizing channels have the greatest minimal output 
Renyi entropy S™ 111 among all maps of the same Renyi entropy of a map 

^map 

One can try to use the extremal position of depolarizing channels to 
estimate the minimal output entropy of some channels. In the case of 
Hastings' counterexample for the additivity conjecture the author showed 
that due to a maximally entangled input state one can obtain smaller 
output entropy of the product of two channels than in the case of any 
product state taken as an input. Let us estimate the Renyi q — 2 output 
entropy for a product channel when the input is maximally entangled. 
Following proposition proved in |51| presents one of estimations. 

Proposition 11. For any entropy S which is subadditive the following 
inequality holds 



|«S map (#i)-5 map ($ 2 )| < S(pi®$ 2 ](|0+><0+|)J < S map ($i)+<S map ($ 2 ), 

(208) 

where \4>+)((j)+\ is a maximally entangled state. 

Proof. The proof starts form the Lindblad inequality 75J, which is based 
on the subadditivity of the von Neumann entropy, 

\S(p) -S p))\<S($(p))<S( P ) + S(^,p)), (209) 

where ?($,p) = [id ®$] M</>)($>|J and \<j)) is a purification of p as in (|106[) . 
The entropy of this state, 5(?(<I > , p)), is the exchange entropy which does 
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not depend on the choice of purification [29] . The state ? defined for 
a channel <3? and the maximally mixed state p* = 1/N is equal to the 
normalized dynamical matrix of $ (|25ll . 

?(<E-,p,) = a 9> = l J D $ , (210) 

The entropy of this state defines S map (<!>). Since the map $ is trace 
preserving, the condition Tr 2 cr<i> = -^-1 holds, see (|26|l . Apply Lindblad 
formula (|209l) to the state 

[*i ® *2](|V+)(V+I) = [*i ® id] ([id®$2](|^+)(V+|)) , (211) 

where 1^+) = -7^ StLi K) ® I*) * s ^ ne maximally mixed state which is a 
purification of p.. Expression (|209p applied to this state gives 

|5 map ($ 2 )-5(?($i®id, <7* 2 ))| < 5(($ 1 ®$ 2 )(|V+)(V+|)) (212) 

< S map ($ 2 ) +5(?($i ®id, a §2 )). 

The exchange entropy S (<;(<3>i ®id, o$ 2 )) is the same as Sk($i, Tr 2 o# 2 )) 
since a purification of a<s> 2 is as well the purification of Tr 2 a<s> 2 and the ex- 
change entropy does not depend on a purification. Due to the trace preser- 
vation formula, Tr 2 o* 2 = p», the state Tr 2 <t* 2 ) = ?('l > i, p*) = 
which completes the proof. □ 

Proposition [11] is applicable for any entropy which is subadditive. The 
Renyi entropy of order q — 2 is not subadditive, however, it is a function of 
the Tsallis entropy of order 2 for which the subadditivity holds. Therefore 
Proposition [11] can be used to estimate the output Renyi q — 2 entropy of 
a product channel if the input state is maximally entangled. The following 
inequality corresponds to Renyi q = 2 version of the lower bound in (|208p . 

-log(l-j e - s ° aP(4,l) - e - s ™ ap (* 2 >j) < S 2 <S» * 2 ) (|^+) (V+ 1)) • (213) 

It is possible to find channels <J?i and $ 2 such that the left hand side of 
(12131) is greater than S™ 111 of depolarizing channel A, which has the same 
map entropy as S™ ap ($i <® <3? 2 ). Notice that for any two channels the map 
entropy of their tensor product is characterized by the following result. 

Proposition 12. The Renyi map entropy 5*™ ap is additive with respect 
to tensor product of quantum maps for any parameter q > 0: 

S™ ap ($i ® $ 2 ) = S7 ap ($i) + S 9 map ($ 2 ). (214) 

Proof. The map entropy S™ ap ($) is defined as the entropy of normalized 
dynamical matrix D<s>. The matrix representation of D& 1 ®& 2 is related to 
superoperator matrix of the quantum operation $i ® $ 2 , due to formula 
(|28[) . Using explicit calculations on matrix elements one can show that 
D<s, 1 ^<s> 2 is unitarily equivalent with D<s> 1 Cg> ZXj 2 . That implies the addi- 
tivity of the map entropies, since the quantum Renyi entropy of any order 
of a given state is a function of its spectrum. 
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Consider a set of iV-dimensional matrices equipped with the Hilbert- 
Schmidt inner product 

(A\B}h := TrA^B. (215) 

In this space the matrix units | i,j= 1,2,..., N} form an or- 

thonormal basis. The elements of this basis are denoted by \i){j\ := \ij)h- 
A quantum operation $ is represented by a matrix $: 

{ij\$\kt) h = Tr(\j)(i\<S>(\k}{£\)), (216) 

hence 

H\k){e\) =J^{i0\H) h \i}{j\. (217) 

*>i 

Due to the reshuffling procedure (|28[l . the entries of the dynamical matrix 
D$ read 

(ab\Dg,\cd) h = (ac\$\bd) h . (218) 
The entries of D$ 1( g)$ 2 are obtained by using unnormalized maximally 
entangled state := Yli e 1*^) ® 1*^) according to definition (I25|) as 

follows 

(obcd|i?* lg |* 2 |e/5ft) = (abcd|[($i®$ 2 )®id](|# + )(*+|)|e/fffc) 

= ^ (aM[(<E > i®4>2)(|^)0'm|)®i^)(jm|]|e/( 7 /i). 

(219) 

Now expression (|217p is used and the matrix elements of D$ 1( g,j> 2 read 
(a6cd|Z?$ 1(g ,0 2 |e/ff/i) = ^ (a^|$i|ij')h (7<5|$2|ij)h (a6cd|«7^> (P5jm\efgh). 

(220) 

Since (afecdlcryi^) is expressed in terms of Kronecker deltas 5 a a5b~,8 C i8di 
and {pSjm\efgh) analogously, the summation over the Greek indexes 
gives, 

(a&crf|.D$ 1 $$ a \efgh) = {ac\D^, 1 \eg) (bd\D<s> 2 \fh) 

= {acbd\D$ 1 ® D* a \egfh). (221) 

The matrix D$ 1( g}$ 2 is related to D^, 1 (g> D<s> 2 by a unitary matrix 
U — Y2 a bed \ a bcd){acbd\. Therefore both matrices have the same eigen- 
values and the same entropies. □ 

Since minimal output entropy of a depolarizing channel A is a function 
of its map entropy (|207|) , the estimation on the left hand side of (|213[) can 
be made in terms of such A for which S™ ap (A) = S™ ap ($i)+S'™ ap (<!>2). As 
a result of this estimation one obtains condition on the pair of channels, 
for which a maximally mixed input state does not decrease the output 
entropy below the smallest value obtained by the product input state, 

MN + 1 g m»P (0l) e -s 2 ro "(0 3 )i < -sf a i > (*ii8* a ) = -[s^Pcso+s^CSa)] 
MN 1 1 ~ 

(222) 
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Figure 8: Colored parts of the figure denote the region described by inequality 
(|222p . This region contains pairs of maps characterized by their map entropy 
for which the additivity is conjectured. The dotted line contains the pairs of 
complementary channels. The region is enlarged if a larger dimensions are 
considered. 




Figure 9: The set of the Pauli channels projected on the plane spanned by the 
map entropy 5' map and the minimal output entropy S*" 1111 . The von Neumann 
entropies are considered. Solid curves correspond to the edges of the asymmetric 
tetrahedron K. The curve AD characterizing the family of depolarizing channels 
does not belong to the boundary of the set. 
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where $1 acts on Mn and $1 on Mm- Fig. [8] presents the region 
defined by (|222p . Such a set is not empty and contains maps, for which 
S™ aP ($i) « S2 nap (*2) or S™ ap ($ 2 ) « 5T P (*i)- The dotted line 
represents the set of complementary channels for which both map en- 
tropies are equal. This set contains the channels breaking the conjecture 
of additivity of minimal output entropy according to the proof of Hast- 
ings. The region defined by (|222|) does not intersect the set. It was also 
shown |40| . |89| that additivity holds if one of the channels is unitary 
or if one of the channels is very noisy. These both cases are covered by 
condition (|222f) . These examples support formulation of 

Conjecture 2 ( |51j). The additivity of minimal output Renyi q — 2 
entropy holds for pair of channels satisfying inequality 1222\ ). 

Recent literature does not answer the question, whether the additivity 
conjecture is broken for low dimensional channels and the Renyi entropy of 
order q — 2. Our Conjecture[3]suggests for which pairs of channels finding 
a counterexample of additivity is unlikely. Conjecture [2] uses the map 



<->l/2 

log(2) 
0.6 

0.4 

0.2 
A_ 



a) 



^5 

log(2)- 



B 



^map 
_Pl/2 



0.6- 

0.4 

0.2- 
A 




2lng(2) 



lo S (2) 



2Iog(2) 



Figure 10: As in Fig. [9h: projection of the set of Pauli channels onto the plane 
spanned by the Renyi entropy of a map 5™ ap and the minimal Renyi output 
entropy S™ ln obtained for a) q — 1/2 and b) q = 5. 



entropy and is formulated for the Renyi entropy of order 2, for which the 
theorem about extremal position of the depolarizing channels was proved. 
This is the key theorem which allows us to derive estimations ()213[) and 
(|222l) . Numerical tests (Fig. [9] and I lOf) suggest that the depolarizing 
channels are not situated at the boundary of the set of all channels in 
the plane (<S'™ ap , S 1 ™'") for q < 2, while their extremal position could be 
confirmed in case q > 2. Nevertheless, the Renyi entropy is a smooth 
function of q. Therefore, a conjecture similar to Conjecture [2] may hold 
also for other values of the Renyi parameter q. 



4.3 Transformations preserving minimal output 
entropy 

In previous chapter the set of one-qubit quantum operations was consid- 
ered in context of the plot (5'™ ap , S™" 1 ). One could ask, whether the family 
of maps lying at the same vertical or horizontal line can be characterized. 
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The following section gives a partial answer to this question. Transfor- 
mations of one qubit maps which preserve the minimal output entropy 
will be considered. Such a transformation changes the quantum chan- 
nel and moves the corresponding point in the plane (S'™ ap ,5'™ m ) along a 
given horizontal line. In the following section we consider the geometrical 
picture of one-qubit maps acting on the set of pure states. One-qubit 
quantum operation transforms the Bloch ball into an ellipsoid inside the 
ball. A transformation of quantum operation which changes the lengths 
of the axes of the ellipsoid and their orientation and leaves the minimal 
output entropy unchanged will be studied. 

Consider the superoperator matrix of a one-qubit quantum operation: 



$ = 



/ 


$11 


$12 


$12 


$14 \ 




$21 


$22 


$32 


$24 




$21 


$32 


$22 


$34 


Vi 


- $11 


-$12 


-$12 


1 - $14/ 



(223) 



Parameters $n and $14 are real, the complex conjugation of $ij is de- 
noted by $ij. The form (|223p guarantees that the dynamical matrix of $ 
is Hermitian and the trace preserving condition (|26p is satisfied. 

Assume that the quantum operation $1 has the output entropy mini- 
mizer at the point 



Pp 



p Vp^-p) 

y/p(l -p) 1-p 



(224) 



Such an assumption is not restrictive since one can always treat the oper- 
ation $1 as a concatenation of a given operation with a unitary rotation 
which does not change the minimal output entropy. The quantum opera- 
tion (|223[) acting on a pure state 
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1 - a 



(225) 



gives an output state 
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which attains the minimum entropy if a — p. 

• Transformation changing the lengths of the axes of the ellipsoid. 

Consider a quantum operation & e iup S oid, which transforms the Bloch 
ball into such an ellipsoid that the end of its longest axis touches the 
Bloch sphere in the "North Pole", 



"ellipsoid 
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V1+V2 

2l|22 
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'71 ~V2 
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1-IJ3 \ 




»?3 



(227) 



• r >9 



Suitable rotations of the Bloch ball before and after the action of 
^ellipsoid guarantees that the point of contact with the Bloch sphere 
is the minimizer of $1. Therefore the concatenation of $1 • ^rotation ■ 
^ellipsoid ■ ^rotation has the same minimal output entropy and the 
same minimizer that $1. The rotation operation is given by 



/ P -sJ{^-p)P -V( 1 ~P)P 1 ~P 

y /(l ~ P)P P P-l -\ /Q--p)p~ 

Vi 1 - p)p p- 1 p -VO- -p)p 

V i-p -p)p Vi 1 -p)p p 

(228) 

where p is defined by the minimizer of output entropy for <I>i . This 
transformation changes the lengths of axes of the ellipsoid but it 
does not change the point at the ellipsoid which is the closest to the 
Bloch sphere. In other words, this transformation does not change 
the directions of the axes of the image of <E>i into the Bloch ball, but 
only their lengths. 

Transformation changing directions of the axis. 

The next transformation changes directions of axes of an ellipsoid 
but preserves the entropy minimizer. In particular, if the image 
of the minimizer is on the longest axis of an ellipsoid, after the 
transformation the point which is the closest to the Bloch sphere is 
no longer on the main axis of the ellipsoid. 

Entropy of an output state (|226[> is a function of its determinant. 
The minimum of the determinant determines the minimum of the 
entropy. Consider a transformation which preserves the value of the 
determinant and compute its derivative in a minimizer. It is useful 
to introduce the compact notation of Eq. (|226[) : 



Pout = aA + (1 - a)B + y/a(l - a)C, (229) 



where matrices A,B and C correspond to the matrices (|226[) . Con- 
sider a transformation $i — » 4>i + ^direction- The output of $i + 
^direction is given by 

d = a(4+i^Sx) + (l-a)(B+i-^=X)+^(T^)(C-X), 
2 ^JP 2y/l-p 

(230) 

where X is a matrix, which is hermitian and has trace equal to 
zero. Moreover, the matrix X satisfies the condition guaranteeing 
that $i + ^direction is completely positive. The state p' coincides 
with (|229[) if a — p. Moreover, the derivative of formula (|229[) with 
respect to a is the same as the derivative of Eq. (|230[) at the point 
a — p. Therefore, the determinants of (|229|) and (|230|) are the same 
and the derivative at a — p is equal to zero. A proper choice of 
parameters in X guarantees that there is a minimum at point a — p. 
Hence both maps, $i and <3>i + ^direction have the same minimal 
output entropy. 
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The part ^direction can be characterized by two parameters (t,n), 




(231) 

Such a form guarantees that the output state of $1 + ^direction is given 
by Eq. (230)) . 

The map $2 of the same minimal output entropy as $1 obtained by 
joint action of three transformations, ^rotation, ^ellipsoid and ^direction, 
on $1 can be given by: 

$2 = <&l<&rotation ' ^ellipsoid ' ^rotation + ^direction- (232) 

We are not able to prove that this transformation contains all possibilities 
of obtaining maps with the same minimal output entropy as a given one, 
however, the transformation is characterized by 5 parameters and also 5 
parameters are needed to have all different (up to one rotation) ellipsoids 
tangent to the sphere on its inner side in a given point. Three parameters 
are associated with the lengths of axes \r}\\, \rfi\, (r/s | , while two parameters 
define the direction of the longest axis n,t. 

Above considerations introduce a 5-parameter transformation of a 
quantum map $1 — > $2- The transformation preserves the minimal out- 
put entropy. Therefore, it determines the family of maps which are situ- 
ated at the same horizontal line of the plot (S" 1 ^ , S min ). Characterization 
of the family of quantum maps parametrized by the minimal output en- 
tropy can be useful to further investigations of relations between <> mm and 
S" 11315 and their consequences. 

5 Davies maps for qubits and qutrits 

Explicit description of general continuous dynamics of an open quantum 
system is difficult in practice. Exact formulas describing the time evolu- 
tion are known in some special cases only. One of the cases in which the 
problem can be solved uses the assumption of a week coupling |91| of a 
low dimensional quantum system interacting with much bigger reservoir 
in the thermal equilibrium. Such an interaction changes only the state of 
the system whereas the state of the environment remains unchanged. By 
analogy to the classical process, in which the evolution of a state does not 
depend on the history, such an evolution is called a Markov process. 

However, while analysing the continuous evolution of the input state, 
sometimes there is no need to know the entire time evolution since only 
the output state is relevant. The "black box" description is useful in 
such cases. A "black box" acts like an evolution discrete in time and can 
be described using completely positive maps, represented as matrices of 
superoperators. 
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Figure 11: Model of a quantum TV-level system characterized by Hamilto- 
nian H interacting with a much larger environment in a thermal equilibrium at 
temperature T. 



The following chapter distinguishes a concrete class of physical pro- 
cesses described by a Davies map [92] . Such a process is compatible with 
the interaction of a quantum state with an environment in a given tem- 
perature, see Fig. 1111 Due to a suitable choice of the entries of a su- 
peroperator matrix $ and relations between them one can say whether 
some continuous time evolution is described by a given discrete quantum 
map. The solution concerns the maps acting on one-qubit, N = 2, and 
one-qutrit, N = 3. In the case of one-qubit maps we determine the state 
which is the most resistant on Davies channels. It will be shown that the 
maximal output 2— norm of Davies maps is additive with respect to the 
tensor product of two such maps. 

5.1 Quantum Markov process 

The quantum Markov process is characterized by quantum maps belong- 
ing to the one-parameter completely positive semigroup, 3> t = expQt, 
where Q denotes a generator and positive parameter t is associated with 
time. 

The most general form of the generator of a completely positive semi- 
group was given by Gorini, Kossakowski, Sudarshan [93] and Lindblad 
[43] . It can be written as 

Q = i8 + C, (233) 

where 5, given by the commutator with the effective Hamiltonian of the 
system S : p — > [p,H], describes the unitary part of the evolution. The 
dissipative part C has the Lindblad form 

C:p^^2(K a pK^ -^{K^K a ,p}), (234) 

a 

where {A, B} — AB + BA is anticommutator, while operators K a can be 
associated with the Kraus representation of the quantum operation. 

Deciding whether a given superoperator matrix belongs to the com- 
pletely positive semigroup was shown [94] to be a problem 'NP' hard with 
respect to the dimension N. Nevertheless, some additional assumptions 
allow one to characterize matrices from completely positive semigroups at 
least for a few low dimensions. In following chapter, such a solution will be 
given for N — 2, and N = 3, under additional conditions: independence 
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of unitary and dissipative parts of the evolution and the detailed balance 
condition. These three conditions define the so-called Davies maps [92 . 
Sometimes the uniqueness of the invariant state is also added to the defi- 
nition. 

5.2 Characterization of the model 

Consider a quantum N - level system characterized by the Hamiltonian 
in its eigenbasis, 

JV 

H = Y,«m\- (235) 

i=l 

Assume that such a system is weekly coupled to the environment of a 
given temperature T, see Fig. [TT] An interaction with the environment 
preserves one invariant state, which is the Gibbs state 

pp = ^exp(-/3H), (236) 

where Z = J^- =1 exp(— /3ei) is a partition function and & — Here 
k represents the Boltzmann constant. A quantum map $ satisfies the 
detailed balance condition if it is Hermitian with respect to the scalar 
product defined by the Gibbs state 

Tr ppA$*(B) = Typd^'{A)B, (237) 

where A and B are arbitrary observables and $* the quantum operation 
in the Heisenberg picture. Detailed description of this condition can be 
found in |96j . 

The name "detailed balance" was taken from the theory of stochastic 
processes. Detailed balance means that in an equilibrium state any two 
levels of the evolving system remain in an equilibrium: the rate of tran- 
sition from the level i to j and the transition rate from j to i are equal. 
Mathematical formula describing this fact reads 

FijPi = FjiPi, (238) 

where Tij are entries of a stochastic transition matrix and pi represent 
the components of the invariant probability vector. 

5.3 Matrix representation of Davies maps 

One qubit map in the "black box" description is represented by a su- 
peroperator matrix. It is a matrix acting on the vector formed by the 
entries of a density matrix ordered in a single column. A superoperator 
$ represents a Davies map, if the following conditions are satisfied. 

• The map $ is completely positive. 

This condition is guaranteed if the Choi-Jamiolkowski matrix D$ 
(|25[) of the map is positive. One has to reshuffle the elements of 
the matrix $ according to ()28|) and check positivity of the resulting 
dynamical matrix = <f> R . 
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Superoperator $ belongs to the semigroup of completely positive 
maps. 

This is equivalent to existence of a generator Q of the Lindblad form 
(|234[) and the parameter t > such that Qt — log$. Knowing the 
logarithm of $ one has to determine whether it is of the Lindblad 
form. It was shown in [95] that if the Choi-Jamiolkowski matrix 
of a given generator is positive in the subspace orthogonal to the 
maximally entangled state, then the generator can be written in the 
Lindblad form. 

It is not a trivial task to write an analytical expression for the loga- 
rithm of a given matrix if its dimension is greater than two. Such a 
problem for 3x3 stochastic matrices is discussed in the last section 
of the following chapter. 

Since the rotational part of the evolution is independent of the dis- 
sipative (contractive) part, the structure of the superoperator is re- 
stricted to the block diagonal form. Off-diagonal elements of the 
density matrix are just multiplied by numbers, while the diagonal 
elements can be mixed between themselves. More detailed discussion 
on this property is given in Section f5. 71 

The detailed balance condition introduces further restrictions on the 
elements of the block acting on the diagonal part of the density 
matrix. This block is a stochastic matrix, the entries of which satisfy 
Eq. p38)) . 



Since now, only the dissipative part of the evolution will be considered. 
Due to the above conditions the dissipative part of the generator of the 
one-qubit Davies maps can be written as 



a,X,p — 



(~a a^-\ 

A 

A 

\a -a^J 



(239) 



while the corresponding superoperator acting on two-dimensional states 
(in the Hamiltonian basis) has the form 



fl-a a^- p \ 

c 

0c 

\ a l-a^-J 



Here, p is a function of temperature, p = (1 + exp (- 
determines the invariant state 



$a, c , P (p*) = p* 







(240) 

(W))" 1 ' which 
(241) 



Notice that Q240p has a block diagonal form which is a consequence of 
independence of rotational and contractive evolution. This is also equiv- 
alent to independence of changes in diagonal and off-diagonal entries of 
a density matrix. The detailed balance condition (|238p implies the form 
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of the outer block in Eq. (|240|) . One-qubit Davies maps form a three- 
parameter family characterized by (a,c,p), where p is a function of the 
temperature. Conditions that such a matrix is an element of the semi- 
group of completely positive maps introduce the following restrictions on 
the parameters (a,c,p): 



a+p<l, 0<c<Jl-- . (242) 

V 1 ~P 

Equality <F = exp£t allows one to write explicit formulas for time 
dependence of parameters a and c, 

a= (1 -p)(l-exp(-At)), c = exp(-n), (243) 



where A and F are parameters such that A > |r > 0. The entire paths 
of the semigroup are showed in Fig. [12] 




Figure 12: Panel a) contains the region of parameters (a, c) allowed by rela- 
tion (|242p and describing the one-qubit Davies maps. The upper border lines 
are also drawn for different temperature T. Panel b) shows the region allowed 
for a given temperature T. The lines describe two semigroup corresponding to 
two different randomly chosen generators C± and £2- The extremal lines corre- 
sponding to the solid line on panel a describe the semigroup with the smallest 
ratio of decoherence to the damping rate. Maps id, $0, ^cg are the identity 
channel, completely depolarizing and coarse graining channel respectively. 

One-qubit Davies map can be written using the Bloch parametrization 

0a 



/ 1 



\ K 3 





m 








'/1 







w 



(244) 



where \r]i\ denote the lengths of axes of the ellipsoid and ii is the translation 
vector. These parameters are related to the parameters (a,c,p) 



7/1 = c > 0, r) 3 

Kl = K2 = 0, 



- -, — >o, 
1-p 

K3 = a fci>o. 
1-p 



(245) 
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The image of the set of pure states under an action of one-qubit Davies 
map is shown in Fig. 1131 The image of the Bloch ball forms an ellipsoid 
with rotational symmetry. Fig. [T3] presents the image of an exemplary 
one-qubit Davies map for which 7/1 > 773, however conditions (|242l) admits 
also the case 773 > <q\. 




Figure 13: Ellipsoid obtained by an action of a one-qubit Davies channel on the 
Bloch sphere. The channel is characterized by parameters (771, 773, K3) defined in 
(12451) . 



5.4 Physical examples 

Qubit maps of the structure similar to (|244[> were analysed before in con- 
text of quantum optics. The unitary evolution is induced by the laser 
field, while the dissipative dynamics is caused by an interaction with the 
environment. The state of a two level atom is characterized by the Bloch 
vector (x, y, z), where z represents the difference between the diagonal en- 
tries of a density matrix equal to the inversion of populations of the atomic 
levels. Variables x and y are associated with the atomic dipole operators. 
The evolution in this set has been defined by means of variables describing 
the decay rate r\ of the coherences and the rate T2 of attaining the equi- 
librium state. These parameters correspond to the variables considered in 
the Section T5. 31 771 = exp(— t/ri) and 773 = exp (— t/rz) which are related 
to squeezing of the axes of the ellipsoid. Formula (|242|) corresponds to 
the relation between the decay rates: 

n < 2r 3 . (246) 

This relation was obtained by analysing a concrete physical model of the 
evolution of the two level system by means of Bloch equations [97| . The 
one-qubit operations (|244l> were also studied by |98j . 
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5.5 Minimal output entropy of Davies maps 



In context of transmission of quantum information, it is natural to ask, 
which pure states are the most resistant with respect to the changes caused 
by the Davies maps. The answer depends on the selected measure of 
decoherence. Such a measure can be described, for example, by means of 
some matrix norm of the output state maximized over the input states. 
Among quantities measuring the decoherence, the minimal output entropy 
is of special importance because some questions concerning the channel 
capacity, such as additivity problem, can be related with similar problem 
written in terms of minimal output entropy. The minimal output entropy 
is related to the maximal norm of the output state if the input is pure. 

Since a Davies map has rotational symmetry, the minimizer can be 
chosen to be a real state: 



H v 
v 1 — /Lt 



(247) 



where v 2 = (1 — fi)fi since the state is pure. After an action of the 
operation (|240[) this state is transformed into 



, _ ((l-a)fj, + b(l- n) cv \ 

p -\ cv a/x+(l-&)(l- M )J' (248) 

where b = ap/(l — p). Computing the eigenvalues and minimizing the 
entropy over fi one can characterize the minimizer in two cases: 

• If c 2 < (1 — a ~ 6)(1 — 26) the minimizer is characterized by fi = 
and it forms an eigenstate of the Hamiltonian H. 

• K c 2 > (1 - a - 6)(1 - 26) the minimizer is characterized by 

(q + 6-l)(2b-l)-c 2 
M 2(a + 6-l) 2 -2c 2 ' 1 ' 

It is no longer the eigenvalue of the Hamiltonian, however, after some 
time of the evolution t ^> the second case changes into the first 
one and the minimizer is a state diag(0, 1). This is an eigenstate of 
the Hamiltonian. The situation that the minimizer is in the vector 
diag(0, 1) reminds the classical evolution of two-dimensional vector 
governed by the stochastic matrix. In this case the extremal vector 
like (0, 1) is the minimizer of the Shannon entropy of the output. 



5.6 Multiplicativity of maximal output norm of 
one— qubit Davies maps 

As discussed in the introduction to Chapter[4] the question of additivity of 
minimal output von Neumann entropy with respect to the tensor product 
of quantum operations is one of the most interesting problem in quantum 
information theory. This problem can be equivalently stated in terms 
of channel capacity. In general, the conjecture on additivity of channel 
capacity is false, however, there is still an interesting problem, for which 
class of maps the conjecture can be confirmed. 
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Recent studies of the problem use the notion of the Renyi entropy 
of order q. This entropy tends to the von Neumann version as q — > 1. 
The problem of additivity of minimal output Renyi q entropy is directly 
related to multiplicativity of the maximal output Schatten q-norm. This 
norm is defined as 

||X|| 9 s = (Tr|Xn 1/9 , (250) 



where \X\ — yXXK. Maximal Schatten q norm of a quantum map $ is: 

||$||f :=max(Tr|$(p)| 9 ) 1/9 , (251) 
p 

where maximization is taken over the entire set of density matrices p. The 
Renyi entropy of order q of a state p can be defined as follows [107] 

■^(pHy^loglHlf. (252) 

Due to logarithm in this formula the multiplicativity of maximal g-norm 
is equivalent to the additivity of minimal output entropy 5*™'". 

In this section, multiplicativity of operator 2-norm induced by the 
Euclidean vector norm will be proved for the quantum one-qubit Davies 
maps. This vector induced norm is not related to the Renyi entropy by 
such an elegant formula like it is in the case for Schatten norm, however, 
it is a bit easier to calculate than the Schatten counterpart. These par- 
ticular results support the general solution for multiplicativity problem 
for Schatten 2-norm which implies the additivity property for minimal 
output Renyi entropy of order 2 and which has been already proved for 
general one-qubit quantum operations |99| (see also |18j). 

5.6.1 Outline of the proof of multiplicativity 

The Euclidean norm (2-norm) of a vector x = (x\, x n ) is defined as: 



£ M 2 . (253) 
\ i=i 



This vector norm induces the 2-norm of an operator A: 

\\A\\ 2 =max^ i^. (254) 



One of the property of this norm (see |100| ) is that ||A||2 is equal to square 
root of the spectral radius of A* A or equivalently to the greatest singular 
value of the matrix A, 

\\A\\ 2 =[r(AU)] 1/2 , (255) 

where a spectral radius r(A*A) — max, and are eigenvalues of A* A. 
In this section the maximal two norm of the output of a quantum map 
$ : Mn — > Mn will be considered 

M* = max ||$(p)|| 2 = max H^ 4 ^ 2 . (256) 
pe-Mjv" "' a>o TtA 
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One can ask, whether the maximal two-norm is multiplicative in a sense: 

M*®n =M 9 M n . (257) 

It will be shown that if <3> is one-qubit Davies map and is an arbitrary 
quantum map acting on ^-dimensional state the multiplicativity holds. 

The idea of the proof of the theorem given below is borrowed from 
the paper of King and Ruskai |54| . These authors prove an analogical 
theorem about a bistochastic quantum map <3>. They noted that the same 
proof holds as well for stochastic one-qubit maps. Here we will present 
an explicit calculations for the case of Davies maps with | rjs | < \t]i\. 

Theorem 8. Let <E> : AI2 — > M2 be an one-qubit Davies map and fl : 
Mn — > Mn be an arbitrary quantum map. The maximal two norm of 
the output is multiplicative: 

Mqqsi = MoAfo. (258) 



In this section the sketch of the proof will be given, while some details 
of the calculation will be presented in the next section. In order to present 
the proof we need to introduce the following set. An arbitrary density 
matrix on H2 ®Hn can be written as a block matrix 




(259) 



where pi, p2, 7 are N xN matrices and the trace condition Tr(pi+P2) = 1 
is satisfied. The output state of the product of two quantum operations 
$ (g> f2, can be described by: 




(260) 



Here 3> denotes an one-qubit operation, while the map Q acts on Mn- 
Also other block matrices will occur and their positivity will play an im- 
portant role during the proof. The Schur complement lemma [101] ensures 
positivity of block matrices, see Lemma [1] Section [3]T] 

To demonstrate additivity (|258|) we shall analyse the inequality M$®n > 
M$Mn which is almost immediate since the equality is attained by a prod- 
uct of states which maximize output norm of each map. Because the entire 
set of states is larger, it contains product and entangled states, the result 
of maximizing over the entire set can give only a better result. Therefore 
to prove multiplicativity of maximal output 2-norm with respect to the 
tensor product of two maps it is enough to show that 

z > M$Mn =^ zl- ($® > 0. (261) 

Insert the block matrix form (|260l) to (|261[) . Due to the Schur complement 
lemma the right hand side of (|261[) is positive if and only if 

Lizl-Py 1 ^ < zl-Q. (262) 
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Notice that this inequality holds if 



\\LL^\ 2 <(z-\\P\\ 2 )(z-\\Q\\ 2 ), (263) 

since using the general property P < \\P\\ one gets: 

Lizl-Py 1 ^ < L{z- ||P|| 2 ) _1 L t < (z-||P|| 3 )~f264) 

< {z- IIQHa) < zl-Q. (265) 

Therefore the positivity of (z — ||P|| 2 ) an d (z — IIQII2) an d inequality (|263|) 
are the only relations needed to prove Theorem [8] These relations will be 
proved in the next section for the case of $ being an arbitrary one-qubit 
Davies map with |r^| < 

5.6.2 Details of the proof of multiplicativity 

Proof, of Theorem [8] It is necessary to find the specific form of M$ , P 
an Q in [[21)0)1 , then to check positivity of {M<s>Mq - \\P\\ 2 ) and (M&Mn - 
HQ || 2 ) j an d finally to prove (|263[l . Let us restrict our considerations to the 
case of Davies maps $, for which 773 < r/ 2 in (12441) as discussed in Section 

Mm 

• Maximal 2-norm of the output, M<s>. 

Use the Bloch parametrization of $ as in (12441) . let it act on the Bloch 
vector (1, x, y, z)\ where x, y, z are real. Moreover x 2 + y 2 + z' 2 = 1 
guarantees restriction to pure states. It is enough to take pure input 
state because the 2-norm is convex on the set of density matrices 
and it attains maximum at the boundary of the set. The spectral 
radius of the square of the output state reads according to Q255p : 

VHHpVHp))) = \i 1 + \/(«3 + zm) 2 + (1 - z 2 H)- (266) 

Since the image of the Davies map has rotational symmetry, there 
are no parameters x and y in this formula. Second derivative of 
the function (|266[) with respect to 2 is negative under the condition: 
V3 ^ r ii- Therefore function (|266|) has a maximum: 

M $ 4(l + y / ^Ml). (267) 

• Output of a product map. 

Now the explicit form of matrices P, Q and L of the output state 
(|260[) will be given. Consider an one-qubit input state. A vector 
(1, x, y, zy corresponds to the density matrix: 

p=\ ( 1 + Z X + iv ). (268) 
2 li- ly 1 — z I 
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Its image with respect to a Davies map (|244p reads: 



Hp) = l + + vi(x + iy) \ (269) 

2 V rji \ x - W l-zrj3-K 3 J 

In the analogous way the initial state in a space Min can be given 
by (259} 

p= l ( * + + * i + A (270) 
2 \ x - ty pi + p 2 - z J 

where z — p 1 — p 2 and x — iy — 27 are N x N matrices. The output 
state of a map <£> ® 1 is: 

(*®id)(p) = (271) 

|(pi +P2 + Vi(pi - P2) + ki(pi + p 2 )) mi 

Vll^ + P2 - %(Pl - P2) ~ K 3 (pi + p 2 )) 

Finally the matrices P, Q and L are defined by comparison of suit- 
able blocks of two block matrices: 

(*®n)W= (ft q) (272) 

in(pi + P2 + %(pi - P2) + k 3 (pi + P2)) ^1^(7) 

»?l^(7) t + P2 - ?73(Pl - P2) - K3(P1 + P2)) 

Multiplicativity. 

One can use the property ||S7(p)|| 2 < Tr(p)Afn (|256|) to show that 
(Af$Mn — ||P|| 2 ) is positive. It is so if 



^M n (l+Jr)l + -^^ j > lMn[Ti(p 1 +p 2 )+r ]s {p 1 -p2)+K3(pi+i>2) 

(273) 

Notice that Tr(pi -+- p%) = 1. To show that the above inequality is 
true, it is sufficient to prove: 



Vi+ f^ 1 2 >^3 + ^ 3 . (274) 
Taking the square of both sides one gets the expression: 

(kzVs ~ (V 2 i ~ Vl)) 2 > 0. (275) 

This implies that (M$M« — II-PII2) > 0- ^ n a similar way we prove 
the positivity of (M<s>Mq — ||Q|| 2 ). The last step is to prove in- 
equality (|263|) . Consider a positive block matrix (l (g> fi)(p) = 



Assume that f2(pi) > (if fi(pi) > one 



71 



can add el to pi and eventually take the limit e — > 0). Due to the 
inequality fi(pi) < ||S7(pi)||2 one can write 

HO( 7 M 7 )V) < \\n( Pl )h{v\n(jMp 1 )- 1 Q(^\v). (276) 

Due to the Schur complement lemma we have fi(p2) > fi(7)n(pi) _1 fi(7)' 
and therefore, 

l|fi( Pl )||2(«|o(7)«(pi) _1 «(7)V> < \Mpi)h{v\n( P 2)\v) < ||n(pi)|| 3 ||n(p2)||a. 

(277) 

Hence the inequality || £7(-y) 0(-y) + 1 1 2 < ||n(pi)||a||n(p2)||2 holds. This 
inequality together with definition (|256l) implies 

||n(7)n( 7 ) t || 2 < M^TrpiTrp 2 . (278) 

Denote Tr(pi) by a:. To prove inequality (1263 II it is enough to show 
that the second inequality holds in the expression below 

ILL*] = n\ H^WIIa < rftx(l-x)M& < (M*Afn-||P|| a )(Af#M n -||Q|| 2 ), 

(279) 

and this is true if 

t£x{1 - x)M 2 a < -Ml Ul + - (r?3(2x - 1) + k 3 ) 2 ] . (280) 

This inequality can be shown by taking the function which is the 
difference between the right hand side and the left hand side. The 
second derivative of this function is equal to 2(r/f — rf^). Therefore 
whenever (rjf > 7/f) the difference is a convex function which has 
minimum at 0. That finishes the proof of the last inequality. There- 
fore inequality (|263|) holds and it proves Theorem [8] 

□ 

In the case < 1 773 1 the proof goes analogously. The maximal output 
norm (|267l) has in this case a simpler form, since the maximizer is a pure 
state described by the Bloch vector (1 = 0, y — 0, z — 1). The specific 
form of the Davies map was used in this proof in (|267p when the formula of 
the maximal output norm was computed and in formula Q271[l . Moreover, 
positivity of K3 is used in (|274[) . 

5.7 Davies maps acting on qutrits 

In this chapter a characterization of the Davies maps for qutrits, iV = 
3, will be given. Going to higher dimensions demands more abstract 
and systematic approach than in the case of one-qubit maps. The entire 
evolution consists of the unitary part and the dissipative part and such 
is the structure of the generator Q — iS + C. The unitary evolution 8 
is governed by the Hamiltonian which in its eigenbasis has a form H = 
Y^ii=i e »l*)(*li where ei > 62 > £3. Differences of energies {wy = e t — ej} 
are called Bohr frequencies. They are eigenvalues of the unitary part 
of the evolution, 5 given by p — > [H, p] , while the eigenvectors of 8 are 
for i,j = 1,2,3. Assume that the set of Bohr frequencies is not 
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degenerated beside the zero frequency case, lju- The subspace related to 
the zero frequency is 3-dimensional. Since the dissipative part C of the 
evolution commutes with the unitary part, it has the same eigenvectors 
and therefore it does not couple the non-degenerated subspaces. Thus the 
off diagonal entries of a density matrix are not mixed with the diagonal 
ones, if the matrix is written in the eigenbasis of the Hamiltonian. 

Like in the case of one-qubit maps only the dissipative part of the 
evolution will be analysed. An one-qutrit Davies map has a structure 
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V 

where Ti\, J~zi, J~32 and (ii, fj,2, A*3 parametrize the map. The off-diagonal 
elements are related by the detailed balance formula 

TijPj = TjiPi. (282) 



Here pi determine the invariant Gibbs state (|236p . The Choi-Jamiolkowski 
matrix of (|28ip preserves the same structure: 
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1-^32 



V 



The generator and its Choi-Jamiolkowski matrix have also the same struc- 
ture. 

Block of the superoperator $ of the Davies quantum operation which 
is related to zero frequency space is a 3 x 3 stochastic matrix 



^31 / 



(1 — -^31 — ^21 Fl2 .Fl3 \ 

F21 I-F21- T Z2 , (284) 

■7"31 -^32 1 — Tv& — 3~23 ) 

where J-32, J-31, Ti\ > 0. Due to the definition of the quantum detailed 
balance condition (|237[) the Davies map is Hermitian with respect to scalar 
product {X, Y)p := Trp^X^Y and therefore it has a real spectrum. 
Moreover, the spectrum is positive, since there is real logarithm of the 
matrix $ represented the Davies map. The positivity of the zero frequency 
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block implies that 

Tzi + T;n + Ti\ < 1, Orz,) 

3 - 4(^32 + Fz\ + F2l) + 3(^32^31 + -FaiJai + .F21J32) > 0. 1 j 

The question considered in this chapter concerns explicit analytical rela- 
tions for entries of the superoperator (|281[) . which imply that the super- 
operator represents a Davies map. One of the condition for $ is that there 
exists an exponential form 

T = e ht . (286) 

Operator L is the zero frequency part of the contractive part of the gen- 
erator of completely positive Davis semigroup. It is parameterized as 
follows: 




(287) 



This is only the zero frequency block which satisfies the detailed balance 
condition. The entire dissipative part of the generator is represented by 
a 9 x 9 matrix. Its Choi matrix has on diagonal elements -L32, L31, 1/2i- 
Since the Choi state of the generator has to be positive on the subspace 
perpendicular to the maximally entangled state we need to require that 
£32,-^31, L21 > 0. 

In the next section an explicit calculation of the logarithm of a stochas- 
tic matrix of order three (|284|) is presented. 



5.7.1 Logarithm of a stochastic matrix of size three 

To compute analytically the logarithm of a positive matrix (|284|) one may 
relay on the following construction. As matrix T has the eigenvalues 
{1, x + y, x — y}, where 

y = ± y/2 Tr T 2 - (Tr T) 2 + 2 Tr T - 3, 1 ' 

the logarithm has the form 




log(.F) =log U x + y ET 1 , (289) 



where U is a unitary matrix which transforms T into its diagonal form. 
Let us evaluate log(J 7 ) without computing the matrix U explicitly. One 
can write that 



log(.F) = \ [log(x 2 - y 2 )Z 2 + log(|±|)z] , (290) 
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where 



/O \ 

Z = U i oh/" 1 




(291) 



The matrix J 7 can be given in terms of Z: 



T = 1 - + xZ 2 + yZ. 



(292) 



This relation allows one to compute yZ — (J~ — 1) — (x — 1)Z 2 . Formula 
for Z 2 can be calculated by taking the square of this equation and using 
the fact that = Z 2 and that Z 2 {T - 1) = (T - 1)Z 2 = ( T - 1). The 
last formula holds since operator (!F — 1) is defined in the subspace for 
which Z 2 is the identity, 



Therefore the logarithm of the matrix T can be expressed according to 
Eq. (|290[) . By comparing a suitable entries of log(J r ) with the parameters 
of L, one gets the parameter L21 as a function of (.F32, .F31, J~2i), 

£21 = 2/2(1 - a; - y) \og{x -y)-(yi(l-x + y) log(x + y)), (294) 

where x and y are given by Eq. (|288|l and 



y 2 ~ 4 - 2y - Tn - J21 + -F13 - J31 + ^23 - J32 + 2 2 " ^ 296) 



The set of points {^23, ^13, ^12}, which defines the set of symmetric 
bistochastic matrices from the dynamical semigroup, is shown in Fig. [14] 
and Fig. Q~5]denoted by E. This set is inside the set of all bistochastic 3x3 
matrices which is denoted by D. The boundaries of the set are stated by 
the constraints I/21, £31, £32 > 0. 

Expression (|294p allows one to check that the set of stochastic matrices 
belonging to the semigroup of completely positive maps with the detailed 
balance condition is not convex. Consider two exemplary points which 
lie near the border of the cross-section and belong to the set (L21 > 
0): {0.5,0,0} and {0.22744,0.22744,0.04512}. Their convex combination 
does not belong to the set. Therefore the set of Davies map is not convex. 
Fig. (represents the cross-section of the set of bistochastic matrices which 
form the zero frequency part of the Davies map represented in the space 
of parameters ^21,^31,^32- Fig. 115b plots a non-convex cross-section of 
the set E by the plane M. 

In order to obtain a full characterization of the Davies map for qutrits, 
not only its zero frequency par have to be analysed. One needs to take 
into consideration also the complete positivity condition and the condi- 
tion on the semigroup related to the Choi-Jamiolkowski matrices of the 
superoperator and its generator. These conditions allow us to specify the 
matrix entries \ii from Eq. (128111 . 

In this way the full characterization of the Davies channels for one- 
qubit and one-qutrits is provided. 



QF-1) \{T- l)-2(s-l)] 
y 2 — (x — l) 2 



(293) 



2y - Ji a - F21 + -F13 - ^31 + J23 - ^32 + 2 



^23^31 



(295) 
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Figure 14: The set E of 3 x 3 bistochastic matrices T (j284|l . which form the 
zero frequency block of the Davies channel $ (|28ip under condition T — > oc, is 
represented by the vector of the off-diagonal elements / = {J-\2, ^23} ■ The 
set E is inside the set of all bistochastic 3x3 matrices D. The characteristic 
points are denoted by / and the corresponding matrix T . 




Figure 15: a) The set E of 3 x 3 bistochastic matrices T (|284jl . which form the 
zero frequency block of the Davies channel $ (|28ip under condition T — > 00, is 
represented by the vector of the off-diagonal elements {J-12, ^13, ^23}- The set 
E is cut by the plane M : J 7 ^ + J~i3 + J~23 = \- The cross-section is presented 
in Panel b) and shows that the set M is not convex. 
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6 Concluding remarks and open problems 



The aim of this thesis was to investigate quantum channels on different 
levels of generality and using different approaches. For instance, general 
properties of quantum channels were considered in Chapter [2l while some 
particular classes of one-qubit and one-qutrit quantum channels were 
analysed in Part IIIII of these thesis. The Davies maps motivated by a 
specific physical model were studied in Chapter [5] Some useful character- 
istics of a quantum channel are provided by different kinds of entropies. 
Among them we used the minimal output entropy, the entropy of a map, 
the entropy of an environment which takes part in an evolution described 
by a channel. Apart of the standard von Neumann entropy which is the 
quantum counterpart of the Shannon entropy, the quantum Renyi and 
Tsallis entropies were also applied. 

In Part [IT] the universal entropic inequality for an arbitrary ensemble 
of quantum states is proved for the von Neumann entropy. This part 
of the thesis treats a quantum channel as a device preparing a quantum 
ensemble. The Holevo quantity of this ensemble is shown to be bounded 
by the entropy of an environment, used in the preparation process. The 
state of the environment after a quantum operation $ is equivalent to the 
output of the complementary channel $. 

One can define selfcomplementary channels for which &(p) = ft(p) for 
any p. Relation (|38f) between the Kraus operators of $ and the Kraus 
operators of $ is useful to specify selfcomplementary channels. Since the 
coherent information (|102[) of such channels is equal to 0, the same holds 
also for the quantum channel capacity (|103[> . Identification of selfcomple- 
mentary channels, as well as investigation of their properties are worth to 
be studied in future. 

Chapter[3]contains the conjecture which establishes a relation between 
the Holevo quantity, and the matrix of fidelities. This leads to a geometric 
characterization of the states in the ensemble. The bound on the Holevo 
quantity proved in Chapter[2]can be also related to other notions of quan- 
tum information theory, such as the quantum discord |102j . |103] which 
measures the quantum correlations in a two-partite system. 

The study of quantum channels is an important task of the modern 
theory of quantum information. For example, the problem of additivity 
of the channel capacity, or equivalently, additivity of the minimal output 
entropy remains open even for channels acting on a single qubit. Re- 
sults presented in this thesis could be further developed to investigate the 
additivity conjecture for different classes of quantum channels. 

Some results of Chapter|4]concern general properties of quantum chan- 
nels. For instance we proved the additivity of the map entropy (|214p . and 
Theorem [6] establishing the extremal position of depolarizing channels in 
the set of all channels characterized by the Renyi entropies 5 , ™ m ('I > ) and 
gmap^-j These results allow us to pose Conjecture [5] specifying pairs of 
maps for which the additivity of channel capacity may hold. 

In Part IIIII some specific types of channels are investigated. Properties 
of one-qubit channels are analysed in Chapter [3] Some transformations 
on one-qubit quantum channels defined in Section 14.31 lead to new re- 
sults on the characterization of the set of quantum channels in the plane 
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(5 ,min ($), S map ($)). The aim of this analysis is to find some conditions 
that enable one to estimate the minimal output entropy, which is difficult 
to compute, by the entropy of the map easy to calculate. 

The Davies channels, which correspond to a concrete physical model, 
are studied in Chapter [S] Superoperators of the Davies maps are specified 
in the case of one-qubit maps and one-qutrit maps. The question whether 
the channel capacity of the Davies maps is additive is still open, although, 
Davies maps acting on JV-level system compose the set of only d = N 2 — 1 
dimensions, while the set of all quantum operations acting on system of 
the same TV has N 2 (N 2 — 1) dimensions. 

The quantum information theory is a modern field of science which 
creates an environment for new future applications and opens new paths 
for development of technology. Quantum channels, which describe any 
possible evolution of a quantum state, play an important role in possible 
applications. Quantum channels describe decoherence caused by the in- 
teraction with an environment. Knowledge of their properties allows one 
to choose the most efficient quantum protocols for a given purpose. Theo- 
retical investigations uncover new possibilities, new laws and fundamental 
restrictions on processing of quantum information. 

The classical theory of information began with investigations on com- 
munication in a given language through given technological tools. How- 
ever, very fast, the laws of information became treated as fundamental 
properties of nature. Therefore, studies in the field of quantum informa- 
tion are so exciting. 
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Appendix 1 



In Appendix we analyze ensembles of three one-qubit states {pi — 
\<t>i}{<t>i\, p2 = \4>2){4>2\, P3} and provide calculations related to Fig. [4] 
necessary to prove Lemma 2] in Section [3.21 

The Bloch vector characterizing the average states can be given by 

OA = o(0,0, 1). (297) 

The Bloch vector representing the mixed state p3 is parameterized by an 
angle a 

OB = 6(0, sin q, cos a). (298) 

The vector OC is chosen in such a way that the ratio \CA\ : \AB\ is 1 : 2. 
Therefore one has 

0~C = ~ (30~A -0~B). (299) 

The point C is in the center of the interval DE, between two pure states 
\cj>i) and |</>2) characterized by the points D and E. Both vectors OD and 
OE form with vector OC the angle 7 so that 

cos7 = |OC|. (300) 

This is in turn the square root of the fidelity \(4>i\<j>2)\ 2 , because the angle 
7 is half of the angle between two pure states, 

-F23 = cos 2 7. (301) 

The fidelity between two one-qubit states represented by Bloch vectors x 
and y reads 

F=~(l+x-y). (302) 
The scalar product of OB and OD is equal to: 

OB ■ OD — b cos (p + 7 — 7) = b ^ cos (p + 7) cos 7 + sin (p + 7) sin 7J . 

(303) 

Hence 



In^.^.J, (OC-OB) 



F 12 = -(1 + OC- OB + b\ 1 - v 4 VI- OC -OC). (304) 

2' V voc-oc v 

The third fidelity F13 can by obtained using Lemma O For /3 = the 
product of three fidelities used in Lemma|4]is a function fa(a, b, a, ft = 0) 
given by 

fo(a, b,a,P = 0) = F 12 F 13 F23 = (305) 
i (9a 2 - 6b cos aa + b 2 ) (b 2 - 3a cos ab - 2) 2 (306) 
+ (9a 2 b 2 (9a 2 - 6b cos aa + b 2 - 4) sin 2 a) . (307) 
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Appendix 2 



In this appendix we present computations necessary to prove Lemma 
[3] It is convenient to change the basis such that the vector OB (see Fig. 
2| is transformed into 

OB = 6(0,0, 1). (308) 

Denote the angle v := n + 7, where /1 is the angle between OB and OD. 
The vectors OD and OE in the new basis can be obtained by rotating the 
state (0, 0, 1) around the axis x by angles: 

0"D = U x (n) (0,0,1), (309) 
0"E = U x (/j, + 2 7 ) (0,0,1). (310) 

The vectors OG and OJ 7 are obtained by rotating the above vectors 
around the axis U x (v)(0, 0, 1) by angle /3. Such a rotation can be defined 
as an action of a unitary matrix U on vectors (|309p . The unitary matrix 
is given by 

U = U z {~ ) U y {v) U z {j3) Ul{v) Ut{~), (311) 
where the rotation matrices read 

(1 \ / cos a sina^ 

cosa -sina , U y (a) = 10 
sin a cos a j \ — sin a cos a / 

(cos a 
sina 


One can use formula (|302[) to calculate the product of three fidelities 
— * / — * / -* / 

for three considered states, OB , OG and OF as a function of the angle 



/ = F 12 F 13 F 23 (312) 
= — cos 2 7 ((cos + cos(27 + (j,) + 2) 2 — cos 2 /3(cos /1 — cos(27 + /i)) 2 ) . 

The product of three pairwise fidelities attains its minimum at /3 = as 
stated in Lemma [3] 
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